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ABSTRACT 



We revisit the causal backreaction paradigm, in which the need for Dark En- 
ergy is eliminated via the generation of an apparent cosmic acceleration from the 
causal flow of inhomogeneity information coming in towards each observer from 
distant structure-forming regions. A second-generation version of this formalism 
is developed, now incorporating the effects of "recursive nonlinearities" : the pro- 
cess by which metric perturbations already established by some given time will 
subsequently act to slow down all future flows of inhomogeneity information. In 
this new formulation, the long-range effects of causal backreaction are damped, 
substantially weakening its impact for simulated models that were previously 
best-fit cosmologies. Despite this result, we find that causal backreaction can be 
recovered as a replacement for Dark Energy through the adoption of larger values 
for the dimensionless 'strength' of the clustering evolution functions being mod- 
eled - a change justified by the hierarchical nature of clustering and virialization 
in the universe, occurring as it does on multiple cosmic length scales simulta- 
neously. With this and with the addition of one extra model parameter used 
to represent the slowdown of clustering due to astrophysical feedback processes, 
an alternative cosmic concordance can once again be achieved for a matter-only 
universe in which the apparent acceleration is generated entirely by causal backre- 
action effects. The only significant drawback is a new degeneracy which broadens 
our predicted range for the observed jerk parameter j° bs , thus removing what 
had appeared to be a clear signature for distinguishing causal backreaction from 
Cosmological Constant ACDM. Considering the long-term fate of the universe, 
we find that incorporating recursive nonlinearities appears to make the possibility 
of an 'eternal' acceleration due to causal backreaction far less likely; though this 
conclusion does not take into account potential influences due to gravitational 
nonlinearities or the large-scale breakdown of cosmological isotropy, effects not 
easily modeled within this formalism. 
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1. INTRODUCTION: COSMIC CONCORDANCE AND CAUSAL 

BACKREACTION 

One of the key questions in Cosmology today relates to the still-unsolved problem of 
what is causing the observed cosmic acceleration. Primarily indicated by Hubble curves 



constructed from lumin osity distance measurements of Type la supernovae (IPerlmutter et al. 



19991 ; iRiess et al.lll998[ ). this (possibly apparent) acceleration is just one aspect of the struggle 
for a consistent picture of the universe; a picture that would also require an explanation o f 



the gap between the observed clustering matter content of f2 M ~ 0.3 (e.g., iTurnerl l2002a[ ) 
and the value of flrnt = 1 as indicated by Cosmic Microwave Background measurements 
of spatial flatness ( jKomatsu et al.ll201ll ). as well as a solution of the "Age Problem/Crisis" 
for matter-only (i.e., decelerating) cosmol ogies in which the universe appears to be younger 
than some of its oldest constituents (e.g., ITurnerl l2002bl ). along with explanations of other 
important issues. 

The standard approach to solving these problems is to introduce some form of "Dark 
Energy" - the simplest case being the Cosmological Constant, A - which can fill the gap 
via Qde = ^Tot — ~ 0-7, which pos sesses negative pressure in order to achieve cosmic 



acceleration (e.g., iKolb fc Turnerl Il990l ) , and which (for non -A cases) recruits some form 



of internal nonadiabatic pressure (e.g., ICaldwell et al.l Il998h in order to avoid clustering 



as matter does. Thus the introduction of Dark Energy (often using spatially-flat ACDM 
models) has led to a broadly-consistent "Cosmic Concordance" - an e mpirical outline which 
has seemed so far to succeed fairly well (e.g., iKomatsu et al.l l201lh at developing into a 
consistent cosmological picture. 

There are serious aesthetic problems with Dark Energy, however, as is well known; the 
most obvious being the problematical introduction of a completely unknown substance as 
the dominant component of the universe. Beyond that, a static (i.e., Cosmological Con- 
stant) form of Dark Energy suffers from two different fine-tuning problems: one being the 
issue that p\ is ^120 orders of ma gnitude smaller than what would be expected from the 
Planck scale (IKolb fc Turnerl 119901 ); and the other being a "Coincidence Problem" (e.g., 
Arkani-Hamed et al.l 120001 ). questioning why observers today happen to live so near the on- 
set of A-domination, given Pa/pm o?- Moving to a dynamically-evolving Dark Energy 
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(DDE), however, invites other problems, since the self- attractive nature of ne gative-pressure 
subst ances (i.e., dE/dV = —P > 0) means that a DDE may cluster spatially ([Caldwell et al. 
19981 ). possibly ruining it as a "smoothly-distributed" cosmic ingredient. This could poten- 
tially be solve d through the ad-hoc addition of some form of nonadiabatic support pressure 
for the DDE (Hulll998l ). but this is a possibility which we have argued against elsewhere 
(lBochnerll201ll ) on thermodynamically-based cosmological grounds. 



Besides Dark Energy, various other methods have been us ed to attempt to explain the 
observed acceleration, such as employing modified gravity (e.g.. lTroddenll2008l ). or assuming 
the existence of an underdense void centered not too far from our cosmic location (e.g., 



Tomital l200ll ). But to avoid the substantial (and perhaps needless) complications which 
arise when assuming departures from General Relativity, as well as the non-Copernican 
'specialness' implied by a local void, we will instead use the feedback from cosmological 
structure formation itself as a natural trigger for the onset of acceleration - a trigger that 
automatically activates at just the right time for observers to see it, due to the fact that all 
such observers will have been created by that very same structure formation which generates 
the observed cosmic acceleration. 



This approach, known generally as "backreaction" , was used by this author in iBochner 



(120 111 ) (henceforth BBI) to find several clustering models which managed to precisely re- 
produce the apparent acceleration seen in Hubble curves of Type la supernova standard 
candles, while simultaneously driving a number of important cosmological parameters to 
within a close proximity of their Concordance values - including the age of the universe, the 
matter density required for spatial flatness, the present-day deceleration parameter, and the 
angular scale of the Cosmic Microwave Background. 

The ability of our models to achieve these goals, despite t he generally pe ssimistic view 
of backreaction typically held by researchers currently (e.g., ISchwara 120101 ) . was due to 
our adoption of an explicitly causal variety of backreaction, which admits the possibility of 
substantial backreaction from Newtonian-strength perturbations. The standard formalism 
used for computing bac kreaction effects, developed through the exte nsive work of Buchert 



and collaborators (e.g., iBuchert fc Ehlersl Il997t iBuchert et al.ll2000l ). is non-causal in the 
sense that it drops all 'gravitomagnetic' (i.e., velocity-dependent) effects, thus rendering it 
unable to account for metric perturbation information flowing (at the speed of null rays) 
from structures forming in one part of the universe, to observers in another. In similar 
fashion, typical studies of cosmic structure formation are also non-causal in that they use 
the Poisson equation without time derivatives of the perturbation potential, thus computing 
metric perturbations from local matter inhomogeneities only, disregarding all gravitational 
information coming in from elsewhere in space. The result is a mistaken (but widespread) 
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notion that the entire Newtonian backreaction Qn can be expressed mathematically as a 
total divergence, thus ultimately rendering it negligible. But by restoring causality with a 
"causal updating" integral that incorporates perturbations to an observer's metric coming 
from inhomogeneities all the way out to the edge of their observational horizon, we find (BBI) 
that the sum of such 'innumerable' Newtonian-strength perturbations - which increase in 
number as r 2 within a spherical shell at distance r from the observer, more than compensating 
for their 1/r weakening with distance - adds up to a total backreaction effect that is not 
only non-negligible (regardless of the smallness of v 2 /c 2 for most matter flows), but is in 
fact a dominant cosmological effect that is fully capable of reproducing the observed cosmic 
acceleration in a fully 'concordant' manner. 

Despite these successes, a major problem with our model in BBI is its utter simplicity: it 
is clearly a toy model, with the results presented there serving primarily as 'proof-of-principle' 
tests, rather than as precision cosmological predictions. Though the simplifications of the 
model are many, one in particular is serious in its consequences, while fortunately being 
not too difficult to fix: specifically, this is the dropping of what we have termed "recursive 
nonlinearities" . Unrelated to gravitational nonlinearities, or to the nonlinear regime of den- 
sity perturbations in structure formation, recursive nonlinearities embody the fact that the 
integrated propagation time of a null ray carrying perturbation information to an observer 
from a distant virializing structure would itself be affected by all of the other perturbation 
information that has already come in to cross that ray's path from everywhere else, during 
all times prior to arrival. In other words, causal updating is itself slowed by the metric 
perturbation information carried by causal updating, creating an operationally nonlinear 
problem. 

This issue was necessarily neglected in BBI, as that work was devoted to introducing our 
'zeroth-order' approach to causal backreaction. But here we fix this problem, incorporating 
recursive nonlinearities into a new, 'first-order' version of our phenomenological model. We 
will find that this alteration significantly changes our results, causing a profound weakening 
of the backreaction effects generated by a given level of clustering, as well as significantly 
damping the long-term effects of information from old perturbations coming in from extreme 
distances. In order to retain causal backreaction as a viable model for generating the observed 
cosmic acceleration - presuming here that this should indeed be done - it will be necessary to 
re-interpret the meaning of our (inherently empirical) 'clumping evolution functions' to now 
consider the effects of hierarchical clustering on a variety of cosmic scales. Doing this, we will 
show that a successful alternative concordance can once again be achieved, with the right 
amount (and temporal behavior) of acceleration, and with good cosmological parameters. 

This paper will be organized as follows: in Section [2j we will re-introduce our original 
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causal backreaction formalism, and then describe the changes implemented in order to incor- 
porate recursive nonlinearities into the model. In Section [31 we will explore the results of the 
new formalism, and discuss the implications of the model parameters that are now needed to 
achieve good data fits. Furthermore, we will discuss how the damping effects due to recursive 
nonlinearities would alter the key factors that determine the 'ultimate' fate of the universe, 
as was discussed for our original formalism in BBI, given an acceleration driven by causal 
backreaction rather than by some form of Dark Energy. Finally, in Section HJ we conclude 
with a summary of these ideas and results, highlighting the role of causal backreaction as a 
fundamental component of cosmological analysis and modeling. 



2. THE CAUSAL BACKREACTION FORMALISM: OLD METHODS AND 

NEW DEVELOPMENTS 

2.1. The Original Toy Model, and its Approximations and Simplifications 

We recall here that the basic premise of the formalism developed in BBI is to phe- 
nomenologically represent the physical processes of structure formation - complex even at 
the level of Newtonian-strength gravitational perturbations - in a simple and convenient 
way. The physics at work within most clustering masses should be as follows: collapsing 
overdensities stabilize themselves and halt their collapse by concentrating their local vorticity 
(or equivalently, by creating a large local velocity dispersion); this concentrated vorticity or 
velocity disper sion leads to real, extra volume expansion in accordance with the Raychaud- 



huri equation ([Hawking Sz Ellisl Il973l ). representable (in the final state) at great distances 
by the tail of a Newtonian potential perturbation to the background Friedman Robertson- 
Walker (FRW) metric; and this Newtonian tail propagates causally outward into space by 
inducing inward mass flows towards the virialized object from farther and farther distances 
as time passes. The total perturbation at time t for any location in space - which will be 
independent of position, assuming similar structure formation rates everywhere - will then 
be the combined effects of innumerable Newtonian tails of this type, coming in towards the 
observer from the virializing masses (in all directions) which by that observation time have 
entered within the observer's cosmological "clustering horizon". 

As is well known, the expansion evolution (i.e., the Friedmann equation) for some spher- 
ical volume V can be deriv ed - using nonr elativistic Newtonian equations, in fact, for a 



matter-dominated universe ( [Weinberg! |1972[ pp. 474-475) - without reference to anything 
outside of that sphere. In contrast, the Newtonian-level backreaction terms which we utilize 
here are due to perturbation information coming in from structures located predominantly 
outside of V (since one must go to cosmological distances for the effects to add up signifi- 
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cantly). For any condensed structure (at great distance) which provides a gravitational pull 
upon the mass in V, its main perturbative effect is simply to impose an extra (Newtonian) 
perturbation potential upon V as an addition to its original cosmological metric. Our phe- 
nomenological approach, therefore, is one in which we model the inhomogeneity-perturbed 
evolution of V with a metric that contains the individually- Newtonian contributions to the 
perturbation potential within V ("$v(0") fr° m au clumped, virialized structures outside of 
V that have been causally 'seen' within V by time t, superposed on top of the background 
Friedmann expansion of V. 

We will reiterate the mathematical essentials of this formalism below, in Section |2T2| but 
first we must recount the various approximations and simplifications which have gone into 
our analysis, to consider their importance and the feasibility of eliminating them in order to 
develop a greater degree of physical realism in these causal backreaction models. 

First of all, though our backreaction-induced metric perturbations will indeed be time- 
dependent (due to the causal flow of inhomogeneity information), they will be entirely 
spatially- independent. As noted above, we do not seek to achieve an observed accelera- 
tion through the mechanism of a local void; but going even further, our model does not 
explicitly include any spatial variations whatsoever. Rather, the system being modeled is 
what we term a "smoothly-inhomogeneous" universe, in which all perturbation information 
blends together evenly in a way that is essentially independent of cosmic position. 

Now, this simplification is one made out of practical necessity, not physical realism. The 
smoothly-inhomogeneous approximation relies upon an assumption of randomly-distributed 
clustering - which is certainly not true, as large clusters are not independent of each other, 
but preferentially clump near one another and are mutually correlated - and this becomes 
ever less true during the ongoing cosmic evolution, as the universe grows more inhomogeneous 
with time. Furthermore, this simplification relies upon the assumption that the region of 
space responsible for the dominant contributions to causal backreaction within volume V 
will be large enough to contain a cosmologically-representative sample of both clusters and 
voids; but as we will see below in Section [31 adding in recursive nonlinearities (to correct 
another simplification, as described below) greatly reduces the size of the cosmological region 
affecting V from what it was in our original toy model, potentially calling this assumption 
into question. 

A proper accounting of causal backreaction in a realistically inhomogeneous universe 
would require the implementation of a fully spatially- detailed, three-dimen sional cosmic 



structure simulation program - perhaps along the lines of ISpringel et al.l ( 120051 ). for example 
- but with Newtonian-level backreaction effects from causal updating now added in. The 
development of such a 3D simulation model is far beyond the scope of this paper (and beyond 
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the efforts of any individual researcher), but would be a useful mission to be undertaken by 
the cosmological community at large. 

A second simplification is our use of Newtoni an potential ter ms - i.e., the long-distance 
approximation of the Schwarzschild metric (e.g., I Weinberg! Il972l ) - to represent the tail of 
each individual perturbation felt from far away , rather tha n using a long-distance approxi- 
mation of the Kerr metric for spinning masses (]Kerrlll963l ). despite the crucial role of some 
form of vorticity in stabilizing most structures against singular collapse. In this case, how- 
ever, the approximation is a good one. A relatively small amount of vorticity can suffice 
for the self-stabilization of a clumped mass, if it is applied perpetually; and the specific 
angular momentum (i.e., [J/(Mc)]) will be small for any mass not on the verge of being an 



extremal black hole. Furthermore, it is easy to show (e.g.. iFranklin fc Bakerl 120071 ) that the 
highest-order deviations from the Newtonian expression in the diagonal metric components 
will go like [J/ (Mcr)] 2 , thus being entirely negligible at the huge distances relevant for causal 
backreaction. (And the leading-order off-diagonal Kerr perturbation terms, though actually 
proportional to (J/r), will effectively cancel out due to angle- averaging in the smoothly- 
inhomogeneous approximation, as discussed in BBI.) Thus it appears quite safe to ignore 
any Kerr-specific perturbation effects for physically reasonable situations. 

Third, our formalism neglects the purely o bservational effec ts of localized inhomo- 
geneities, such as lensing along beam paths (e.g., iKantowskil 120031 ) for rays from standard 



candles, and similar perturbative effects upon the apparent luminosity distance relationship 
for rays passing through inhomogeneous regions, which in some models represents t he pri- 
mary 'backreaction' effect resulting from structure formation (e.g-. lBiswas et al.ll2007l ). Even 
if such effects by themselves are too small to generate an observed cosmic acceleration, they 
will still alter the output parameters estimated while using any cosmological model (includ- 
ing ours), and thus should be kept track of; and in case our causal backreaction method 
also falls short of providing the full result of an apparent acceleration all by itself, it might 
successfully be combined with these other observational effects upon the light rays to pro- 
duce an apparent acceleration once everything is added together (this point to be discussed 
again in Section l3~TTj) . Combining these purely observational effects with those from causal 
backreaction is therefore an important task, and likely quite a feasible one; though not one 
addressed yet in this current paper. 

The next, most theoretically treacherous approximation is our neglect of nonlinear grav- 
itational effects, a simplification made implicitly by our method of linearly adding together 
the individual metric perturbations contributed by different self-stabilized mass 'clumps' 
in order to produce the total, summed, Newtonian-strength perturbation potential. (Note 
that this is not the full "Newtonian" approximation usually employed, since while we do 
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assume weak-gravity, we do not completely assume 'slow-motion' - in the sense of drop- 
ping all time derivatives of the perturbation potential - as that would neglect the causal 
flow of perturbation information.) Thus our formalism explicitly neglects the nonlinear, 
purely general-relativistic effects that most other researchers primarily focus upon when 
studying "backreaction" . This approximation becomes increasingly bad as the magnitude 
of the summed perturbation potential approaches unity; but since (as will be seen below) 
this potential typically does not grow to values in excess of ~0.5 — 0.6 or so as t — > t for 
most of our best-fitting simulation runs, the approximation is probably good enough for our 
simulations to provide fairly accurate estimations of the cosmic evolution up to now, and of 
our measurable cosmological parameters. (And to the extent that it is not good enough, a 
significant contribution due to nonlinear gravitational terms would likely only help produce 
the desired acceleration even more easily.) Thus it is probably not necessary for us to include 
higher-order gravitational terms in our formalism, in order to achieve a sufficiently reliable 
understanding of the currently-observable universe for our present purpose of pointing the 
way towards an alternative concordance; a fortunate situation, since our model is fundamen- 
tally designed around a linearized-gravity approach, and it may be challenging to find any 
convenient way of modifying it to include nonlinear gravitational effects. On the other hand, 
given the ever-increasing strength of gravitational nonlinearities in the cosmos over time, a 
fully general-relativistic model of causal backreaction (computed using a 3D simulation of 
realistically-distributed inhomogeneities) would almost certainly be necessary for accurately 
predicting the long-term future evolution of the universe. 

Lastly, there is the approximation regarding what we have referred to as "recursive 
nonlinearities". As will be seen from the metric given below in Section 12.21 one of the 
effects of causal backreaction is real extra volume creation. But since causal backreaction 
depends upon the propagation of inhomogeneity information through space, the extra volume 
produced by old information from perturbations will slow down the propagation of all future 
inhomogeneity information (as well as carrying all perturbing masses farther away from all 
observation points), thus feeding back upon the causal backreaction process in such a way 
as to strongly dampen it. Of all of the simplifications and approximations discussed so far 
in this subsection, the neglect of these recursive nonlinearities most likely has the strongest 
impact upon the quantitative predictions emerging from our causal backreaction models. 
Fortunately, however, fixing this problem by adding these recursive nonlinearities into our 
formalism is one of the simpler improvements in physical realism for us to make; and hence, 
this paper focuses upon achieving this fix, and then calculating and interpreting the results 
produced by this 'second-generation' causal backreaction formalism. 
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2.2. The Old Formalism and its Results 

Here we recall the technical details of our original formalism developed in BBI, to set 
the stage for its further development to follow. 

To obtain the Newtonian approximation of a single 'clumped' (i.e., virialized, self- 
stabilized) object of mass M, embedded at the origin (r = 0) in an expandi ng, spatially-flat , 



matter-dominated (MD) universe, one may linearize the McVitt ie solution (IMcyittidll933l ). 
as can be seen from the perturbed FRW expression given in iKaloper et al. I J201(j l. The 



resulting Newtonianly-perturbed FRW cosmology is given by the metric: 
ds 2 « -c 2 [l + (2/c 2 )$(t)]dt 2 + [a M D(^ , (1) 

where = {— GM/[aMD(t)r}}, and omd(^) oc t 2//3 is the unperturbed MD scale factor 
evolution function. 

For our model of a smoothly-inhomogeneous universe, we assume a random distribution 
of clustered masses, being the same essentially everywhere and in every direction. So turning 
the above expression around, we consider the situation for an observer at the origin, whose 
metric is affected by a collection of (roughly identical) discrete masses - for now confined 
to a spherical shell at coordinate distance r', and with total mass M - that are distributed 
randomly in direction. The different directions of the various clumps does not matter for the 
observer's g u metric component; but it does matter for the spatial metric components, due 
to the fact that a spatial displacement from the origin would pick up (direction-dependent) 
factors of cos 2 9 in ds 2 from the different angles of the motion with respect to the individual 
clumped masses, since only the radial projection of a given translation (with respect to a 
particular clump) will 'feel' the perturbation potential from that clump in its contribution 
(all within g rr ) to the interval ds 2 . Averaging over direction in three dimensions, the spatial 
metric terms therefore pick up a factor of (cos 2 9) = (1/3) for the total effect when summing 
over all of the discrete masses; and the total, (gravitationally-)linearly summed and angle- 
averaged metric for this observer at the origin can then be written in 'isotropized' fashion, 
as: 

ds 2 = -c 2 {l - [iWt)/r']} dt 2 + [a MD (t)] 2 {l + (1/3) [iW*)M} W , (2) 

where R Sch (t) = {(2GM/c 2 )/[a MD (t)]}, and \df\ 2 = [dr 2 + r 2 d9 2 + r 2 sin 2 9d<f) 2 ) = \dx\ 2 = 
(dx 2 + dy 2 + dz 2 ). (Note that this factor of 1/3 is not 'fundamental', but is merely the result 
of our approximating the linearized sum of many individual 'Newtonian' solutions, which 
effectively spreads out the total spatial perturbation among all three spatial metric terms, 
rather than confining it solely to usual. Thus the general relativistic expectation 

of equal temporal and spatial potentials - i.e., ip = <p ~ is not really violated here, and no 
actual new physics or modified gravity is implied by it.) We regard the term multiplying 
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[omd(£)] 2 |c^1 2 as a 'true' increase in spatial volume, and the g u term as an 'observational' 
term, slowing down the perceptions of observers at any given t. Thus, even if the spatial term 
by itself is not sufficient to generate a real volumetric acceleration, the spatial and temporal 
perturbations coupled together may indeed be enough to create an "apparent acceleration" 
capable of explaining all relevant cosmological observations. 

The expression in Equation [2] represents the perturbations to the metric due to masses 
at some specific coordinate distance r' - and thus from a specific look-back time t' - as 
seen from some particular observational point. The total metric at that spacetime point 
must be computed via an integration over all possible distances, out to the distance (and 
thus look-back time) at which the universe had been essentially unclustered. Finally, a light 
ray reaching us from its source (e.g., a Type la supernova being used as a standard candle) 
travels to us in a path composed of a collection of such points, where the metric at each point 
must be calculated via its own integration out to its individual "inhomogeneity horizon"; 
and only by calculating the metric at every point in the pathway from the supernova to our 
final location here at r = z = can we figure out the total distance that the light ray has 
been able to travel through the increasingly perturbed metric, given its emission at some 
specific redshift z. 

Finite look-back times imply that when one is feeling the effects of clumps at cosmolog- 
ical distances, what one is really sampling is the clustering as it was at an earlier, retarded 
time. Thus the cosmological evolution at all times is dependent upon the entire history of 
the development of clustering. In BBI, we defined a heuristic "clumping evolution function", 
*&(t), intended to simply represent the fraction of the cosmic matter that has completed 
its relaxation to a self-stabilized, steady state by cosmic time t, as opposed to that mass 
still freely expanding (or collapsing) within a still dynamically-evolving patch of space. The 
function \l/(t) was therefore defined over the range from =0 (perfectly smooth matter), 
to \l/(t) = 1 (everything clumped). This is an extremely simplified way of representing the 
virialization of cosmic structures, and we will in fact be forced to re-evaluate the meaning 
(and the range) of the function \l/(i) later on in this paper. But for now, we just recall that 
the specific functions used for \l/(t) in our numerical simulations are not rigorously derived 
from first principles, but rather are adopted for simplicity, with their functional forms and 
input parameters being motivated by physically reasonable structure formation behaviors, 
in conjunction with certain quantitative cosmological measurements. 

Now, to determine how Hubble curves obtained from standard candle measurements are 
calculated in our formalism, consider a light ray emitted by a supernova at cosmic coordinate 
time t = £sn, which then propagates from the supernova at r = tsn, to us at r = 0, t = t . 
We refer here to the geometry depicted in Figure [TJ 
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For each point P = (r, t) of the trajectory, the metric at that point will be perturbed 
away from the background FRW form by all of the virialized clumps that have entered its 
causal horizon. Consider a sphere of (coordinate) radius a, centered around point P, with 
coordinates (a,t ret ) (where t re t < t is the retarded time), defined such that the information 
about the state of the clumping of matter on that sphere at time t re t will arrive - via causal 
updating, traveling at the speed of null rays - to point P at the precise time t. To compute 
the fully-perturbed metric at P, we must integrate over the clumping effects of all such 
radii a, from a = out to a max , the farthest distance from P out from which clumping 
information can have causally arrived since the clustering of matter had originally begun in 
cosmic history. 

For the remaining calculations in this subsection, note that we will be using our old 
formulation from BBI, without recursive nonlinearities (those changes will be presented later, 
in Subsection 12 .3 j) . In particular, this affects the speed of the propagation of inhomogeneity 
information through coordinate space (as well as some other issues), which for now will be 
calculated with respect to the unperturbed FRW backreaction metric, rather than referencing 
the perturbed metric itself in an explicitly self-consistent manner. 

For an FRW metric with a(t) = ao(t/to) 2 ^ 3 ; the coordinate distance traveled by a null ray 
in the cosmic time span from t\ to t 2 will be a = (c/ao) J^(t/to)~ 2 / 3 dt = [(3c/ao)(to) 2//3 (^2^ 3 — 
t\ )]. Defining a = 3ct = 2c/ H , and with t 2 = t, ii = t iet , we thus have: a = [(t/to) 1 ^ 3 — 
(tret/to) 1 ^ 3 ]. We then turn this into a prescription for computing t ret as a function of t 
(relative to the present time, to) and a, as follows: 

t ret {t,a) =t [{t/t ) 1/3 -a] 3 . (3) 

Similarly, we can determine a max , given some initial time ti^t at which structure formation 
can be reasonably said to have started (i.e., ^(t < ti nit ) = 0): 

«max(t, tlnit) = [(t/t ) l/3 - (t lnit / to) 1 f 3 ] . (4) 

Now, how the metric at P is affected by a spherical shell of material at coordinate radius 
a depends upon the state of clumping there at the appropriate retarded time: ^f[t TCt (t,a)]. 
The total effect is then computed by integrating all shells from a = out to a = a max (t, ti m t); 
but in order to compute the metric perturbation from each shell quantitatively, it is first 
necessary to relate this clumping function to an actual physical density of material. 

As discussed above, for now consider the ^(t) function as representing the dimensionless 
ratio of matter which can appropriately be defined as 'clumped' at a given time, expressed 
as a fraction of the total physical density. Assuming a flat MD cosmology as the initially 
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unperturbed state, the total physical density at all times will merely be an evolved version 
of the unperturbed FRW critical closure density from early (pre-perturbation) times. 

Recalling Equation^ we have the perturbation term [Rs c h(t) / r '] = {{2G M / c 2 ) / [r' cimd (£)]}> 
with a MD (f) = a (t/t ) 2 / 3 = c[18t 2 / Hq} 1 ^ 3 . The value of M to use here is given by the 
clumped matter density at coordinate distance a times the infinitesimal volume element of 
the shell. The clumped matter density at time t, as implied above (and still as considered 
with respect to the unperturbed background metric), will equal (t) p crit (t)]; and the vol- 
ume element in the integrand representing the effects of that shell will be 4:TrR 2 hys dR p hy S = 
47r[a M D(£) a] 2 [a M D^) da]. 

Collecting these terms (and letting if = t Tet (t, a)), the integrand will thus be equal to: 



[Rsch(t)/r'} r , =a ^ ia+da) = {(2G/c 2 ) dM I [a MD (t) a]} (5a) 

= {(2G/c 2 ) [omd(*) a]' 1 p crit (t)] [AiiR 2 phys dR phys ]} (5b) 

= {(8nG/c 2 ) [a MD (t) a]' 1 [ Pcrit (t) [a MD (t)] 3 ] [a 2 da]} (5c) 

= {(SttG/c 2 ) auvit)- 1 [ Pcrit (t ) a 3 ] [ada]} (5d) 

= {(SttG/c 2 ) *(t') [(toA) 2/3 (Scto)- 1 } {[3H 2 /(8nG)] (3ct ) 3 } [aejfifc) 

= {12 *(0 [(t A) 2/3 ] [ada]} , (5f) 



where for simplification we have used H = (2/3)t () 1 and the fact that [p(t)a(t) 3 ] is constant, 
both true for a matter-dominated universe. Note also that only \l/ is evaluated at the retarded 
time, t re t (t, en). The strength of the metric perturbation (at time t for point P) from a 
point-like Newtonian perturbation embedded in the expansion actually depends upon its 
instantaneous physical distance from P at t, as is obvious from $ = {— GM/[omd(^) r ]} = 
[—GM/Rpinyaty)] in Equation [TJ The only "relativistic" piece of propagating information 
which is causally delayed is the state of clumping, ty[t Tet (t, a)], that has just then arrived 
from coordinate distance a to observer P at (r, t). 

The total integrated metric perturbation function due to clumping, I(t), as experienced 
by a null ray passing through point P at (r, t), is thus calculated via the causal updating 
integral: 

/'Q'max (i,iinit) 

/(*) = / {12 #[*„*(*,«)] [(V*) 2/3 ]} a da , (6) 

Jo 

with I(t) implicitly being a function of ti nit (with I(t) = for t < iinit), as well as of to- 
Inserting this result back into the formalism of Equation [21 the final clumping-perturbed 
metric that we will use for all of our subsequent cosmological calculations becomes: 

ds * = _ c 2[! _ m] dt 2 + { [ aMD (t)] 2 [1 + (l/3)/(t)]} \df\ 2 . (7) 
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Given this metric, we can define relationships between all of the unperturbed FRW 
("bare") cosmological parameters, and their corresponding observable ("dressed") parame- 
ters. We must first consider cosmological redshifts, where we have: 

^ {t) = ^(M_ 1 = ( t o/t) 2/3_ li (8) 
ClMD{t) 

and: 



where quantities without superscripts like to = ^o RW ; e ^ c -; wm generically refer to unper- 
turbed FRW parameters; and where values of observational quantities in our backreaction- 
perturbed formalism will be expressly indicated (e.g., tQ bs ). 

Now, in order to calculate observed luminosity distances, we must compute the coordi- 
nate distance r of a supernova going off at coordinate time t, for which a light ray would 
be arriving here (at r = 0) precisely at to- For a null ray and pure inward radial motion we 
have: 



rss'(t) = \r^(t )-r—(t)\= I ^} dt> (10a) 




9tt(t') 



to 



^\lTTjmm*« (10b) 



This coordinate distance function can then be converted into an expression for the 
observed luminosity distance by adapting the expression from the homogeneous FRW case, 
for which = [a tsn (1 + z)], as follows: 

d w (t) = [aoVl + (l/3)/(t )] r F ™(t) [1 + z° b s (t)] (11a) 

(/o/3) ctT r Vr 2/3 1 1 t^L> dt (nb) 



y/l + [I(t)/3] J t LV ' V l + [/(tf)/3]- 

l + (/ /3) Ct r\ u ,,-2/3 I 1-I(t'r 



y^m^^jJ^iTTmk^- (llc) 

where Io = /(to), and t r , t' r are dimensionless time ratios (e.g., t T = t/to), with no change 
to the essential form of J(t) (i.e., J(t) = I(t r ■ t ) =>- I(t r )). 
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All cosmologically-relevant curves, fits, and parameters can now be calculated from 
the metric (Equation [7]), and from this luminosity distance function (Equation [TT]) and 
its derivatives, as investigated in BBI. Important modeled quantities include: the residual 
distance modulus function (with respect to an empty coasting universe), A(m — M) Pert (z° bs ), 
and the quality and probability of its fit, Xfh an d Pfh, to the Type la supernova (SNIa) data; 
the observed (versus unperturbed) Hubble Constant, if° bs (versus if FRW ); the physically- 
measurable age of the universe, t% bs ; the "true" value of the total cosmic matter density, 
f2^ RW , which determines the primordial spatial curvature (i.e., f2^ RW = 1 for flatness in the 
pre-perturbed epoch), which our model must relate to some measured value of the density, 
f2^ bs , for normalization (here and in BBI we use ^M bs = 0-27); the observable values (defined 
for z —7- 0) of the deceleration parameter g° bs , the effective (total) cosmic equation of state 
w® hs , and the jerk parameter j° bs ; and finally, for comparison with complementary data sets 
from much earlier cosmic epochs, we compute the acoustic scale of the Cosmic Microwave 
Background (CMB) acoustic peaks, Z° bs . 

As a technical note, we recall that all evolving quantities in our numerical simulation 
program are calculated as discrete arrays in t (and in z° hs (t)), with a tested pixelization 
that is fine enough for great accuracy in all parameters. Given that the discrete version 
of dr 

,Pert(* s ) must be differentiated (and evaluated specifically for z — > 0) to obtain cos- 
mological parameters, we do so by using the definition of the derivative for each pixel, as 
follows: 

, j(AT-l)//, n ,(AT-1)/,, \ 

U N ' 1 \A N - 1 )'-[ = Q L,Pert »L,Pcrt \ T {i} ) ,* \ 

L<VertJ« " d zQhs L^L,Pert J« " ^Obs (%+1}) _ ^Obs^) > M 

and then obtain the z — > limit from the last, latest-in-time pixel: 

[^Pert](z-*-0, t^t ) = [^Pert]{Wp ix } • (13) 

The cosmological results which we obtain from this procedure appear to be robust, with only 
minor difficulties, as will be discussed below; and while the greatest discretization errors occur 
for j® hs , which requires three differentiations of d^ } p ert (t), virtually all of the results which 
we quote here for j^ bs should be well within 1% in terms of numerical accuracy. 

In order to conduct specific calculations with our formalism, we must design a set of 
physically reasonable clumping evolution functions, *&(t), to serve as convenient proxies for 
the combined effects of the linear density evolution of early-stage clustering, the nonlinear 
regime and virialization for very dense clumps, and the initial development (in many cases 
triggered by collisions) of entirely new clumps with substantial mass. Guided by general 
cosmological considerations and simplicity, in BBI we chose three different classes of time- 
dependent behaviors to examine: ^(i) oc a(t) oc t 2 / 3 , whi ch is proportional to the evolving 



contrast of a density variation, Sp/p, in the linear regime (IKolb fc Turnerlll990l ); $?(t) cx t, a 
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generally sensible choice depending simply upon the amount of time available for clumping; 
and ^(t) a t 2 , an 'accelerating' clumping model which we initially chose as a test case to 
see whether that would possibly help in creating an observed acceleration. This last class 
of models will take on more significance in this paper, however, since our results below with 
recursive nonlinearities will force us to regard \l/(t) as not merely a simple percentage of 
clumped versus unclumped matter, but as a quantity reflecting the details of virialization 
brought to completion (i.e., extremely nonlinear density perturbations) on multiple cosmic 
scales; and appropriately, the density con trast evolution for inh omogeneities in the nonlinear 
regime goes as Sp/p oc a(t) n with n > 3 faolb fc Turner 1990 . p. 322) - that is, 5p/p oc t m 



with m = (2/3)n > 2 for a matter-dominated universe. 

Quantitatively, we defined our three different classes of clumping evolution models as 
follows: 

(14a) 
(14b) 
(14c) 



*md(*) 


EE 




^Init ^2/3 




~~ *Init 






*Un(t) 


EE 




— ^Init \ 




— ^Init 




^Sqr(t) 


EE 




— ^Init \2 







Here ti nit (or equivalently, Zhdt = [(^oAinit) 2 ^ 3 — 1], as per Equation^]) represents the effective 
beginning of clumping, such that < ti nit ) ee for all models; and = "if (to) represents 
the current state of clumping today. Each of these models has two physically meaningful 
parameters to vary (besides the universally-optimizable parameter for all models, Hq): ^q, 
for which we selected values by estimating the fractions of Dark and baryonic matter that 
may likely be clumped by to] and zi nit , for which we selected values by roughly linking the 
beginning of clumping with the process of cosmic reionization. 

In BBI, we evaluated 60 different models, using ^ = (0.78, 0.85, 0.92, 0.96, 1.0) and 
^init = (5, 10, 15,25) for each of ^un(t), *mdW, and ^sqr(i)- We found that these models 
produced curves which behaved very much like ACDM, with the \l/sqr runs looking like flat 
ACDM with n A ~ 0.3 - 0.4, the ^ Lin runs looking like fi A ~ 0.5 - 0.8, and the ^ M d runs 
looking like fi A ~ 0.65 - 0.97. 

As was expected, due to the increase in the number of perturbing clumps with distance 
(as r 2 ), the total perturbative effects were typically dominated by the largest distances (and 
thus the earliest look-back times) out to which one could still see significant inhomogeneities; 
hence clumping evolution functions with a more rapid onset of clustering at early times would 
produce a more powerful causal backreaction effect. Thus the $md models were by far the 
strongest, while the \l/sqr models were all too weak to reproduce the observed acceleration, 
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with the ^Lin(^) models falling intermediate between the two; and in all cases, the choice of 
an earlier zi nit led to a stronger overall backreaction. 

All told, roughly a dozen of these models resulted in good fits to the standard candle 
SNIa data, reproducing the observed cosmic acceleration essentially as well as the best- 
fit ACDM Dark Energy model; and about half of those dozen models produced especially 
good cosmological parameters, as well. The overall picture was that of a successful attempt 
at achieving an alternative cosmic concordance without Dark Energy. Additionally, there 
emerged the bonus of a testable prediction which could distinguish causal backreaction from 
a Cosmological Constant: a result of j° bs ^> 1 for all of our models which produced good 
fits, in contrast to the mandatory value of j° bs = 1 for all flat ACDM models, regardless of 
the value of 

With those results in hand, the next concern was to determine how the effects of recursive 
nonlinearities might modify these outcomes, once incorporated into our numerical models. 
Our initial expectation was that the overall Hubble curves of each modeled cosmology would 
likely not change significantly, except perhaps at the end, as z — > 0; a result that would 
not alter the SNIa fitting results much, but which could have a measurable effect upon 
the cosmological parameters evaluated at t - particularly upon our 'paradigm falsifiability' 
parameter, jQ hs . The actual outcome, however, is that the proper treatment of recursive 
nonlinearities ends up altering the situation to a great degree, requiring us to re-evaluate 
(see Section [3]) the ultimate implications of causal backreaction; although it does still remain 
true that an alternative concordance can be successfully achieved with it. First, however, in 
the following subsection we give the precise details explaining how recursive nonlinearities 
for causal backreaction are technically implemented in our numerical simulations. 

2.3. The New Formalism: Incorporating Recursive Nonlinearities into Causal 

Backreaction 

Though in practice necessitating a complete re-write of our numerical simulation pro- 
gram, the essential features of recursive nonlinearities require just two fundamental changes. 

The first change involves the rate at which clustering information can propagate through 
the inhomogeneity-perturbed universe, as depicted earlier in Figure [TJ The result is that 
t Tet (t, a) and a mBX (t, ti Qit ) are altered from their simple functional forms (as per Equations EE]) 
for a matter-dominated FRW universe, and now depend upon the metric perturbation po- 
tential function, I(t). But as we see from Equation [61 1(f) is itself calculated via an integral 
depending upon t ret (t,a) and a max (t, tinit) - Hence I(t) is now defined in terms of I(t), in a 
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fundamentally recursive fashion. 

Second, the extra spatial volume created by nonzero I(t) will have a dilution effect 
upon the strength of the perturbation induced by any (now farther-away) clumped mass 
upon the metric at our observation point. Examining Equation [5] (specifically |5k), we 
recall that the perturbation term for each spherical shell of clumped matter is given by 
[Rsch{t) /r'} r >=a^(a+da) = {(2G/c 2 ) dM / [cimd(^) £*]}. Now, the differential mass element 
in the shell, dM, does not actually change (and therefore requires no correction factor), 
since the dilution of its mass density is precisely offset by its expanded volume. But what 
does change is the effective distance of those perturbing clumps from the observation point. 
Recalling that the strength of a Newtonian perturbation embedded in an expanding uni- 
verse depends simply upon its instantaneous physical distance, we see that the denominator 
[omd(£) a ] must now be replaced by the term: [J~g^- a] = [a MD (t) a/1 + (l/3)/(t) a], which 
in the end just puts a factor of yl + (l/3)/(t) into the denominator of the integral for lit). 

Therefore, the modified formula (replacing Equation [6]) for calculating the metric per- 
turbation function I(t) with the incorporation of recursive nonlinearities ("RNL"), is given 
as: 

jR NL (t) = {i2 nut a, i)} [(t /tm}a da ^ (i5) 



/o y/1 + (1/3) 

where the term "/" inside the integral on the right-hand side represents the actual function, 
J RNL (t), itself. But since the denominator term is in fact independent of the integration 
variable a, we can remove it from the integral and bring it to the left-hand side of the 
equation. The metric perturbation function can therefore be (numerically) solved for any 
given t as the solution of: 



rRNL 



p Ot. m ax 

(t) y/1 + [J RNL (t)/3] = / {12 *[t ret (t, a, I)} [(toA) 2/3 ]} a da . (16) 

Jo 



To evaluate the expressions t ret (t,a,I) and a m3iX (t, ti n i t , /) in this above integral, we 
utilize (analogously with the discussion preceding Equations I3"]H|) : 

(TT n _ 1 f T V 1 ~ /RNL ( r dT ' 

1 ' ^ j 3 7 rrct ^ + [/RNL (T0/3] (T') 2 /3 ' {U) 

with T = (t/t ) and T ret = (t rct /t ). This expression for a(T,T ret , I) could in theory be 
inverted to produce T ret (T, a, I); and in addition, we have a max = a(T, T Init , I). 

In practice, we perform these recursively-defined 'integrals' by utilizing discrete arrays 
in cosmic coordinate time with some large number of pixels covering the range from Ti n j t to 
T = 1, for which the later pixels are calculated in terms of the earlier pixels. Beginning 
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with the first pixel at Ti mt , we thus have the recipe (with iV P i x pixels, and AT = [(To — 
Ti nit )/(iV pix -l)]): 

cw.1 = I RNL = 0, (18a) 
1 AT 

«max,2 = , / 2 RNL = , (18b) 

and then, for i = {3, iVp ix }: 



h 7~R.NL 

1 V J (i-i) AT 



amax,i = Q!max,(i-1) + (q / = 12/3 } > ( 19a ) 



-[2 

^-rnl _ _ __ {\J/[T(j_ fe )][a maXii - a maX) (i_fc)][a maXi (i + i_ fc ) - amax^j-fc)]} (19b) 

T i fc={l,(i-2)} 



Rnl ./i+[/RNL/3]=Xf NL . (19c) 



The result of this iterative loop is the discrete array {/ RNL }, which serves as 7 RNL (t) for 
all cosmo logical calculations, as described above. As a test, we have verified that removing 
the effects of recursive nonlinearities (essentially adjusting to unity all terms inside the radical 
signs in Equations [T9"r.c) succeeds in reproducing the results of our old formalism to a great 
degree of accuracy; and for our new model simulations with the recursive nonlinearities 
included, we have checked to make sure that the full suite of results with this new formalism 
appears quite sensible and consistent in all cases. 

We have also tested the pixelization for precision of the output results, finding that 
runs with ~1000 pixels are sufficient for quick model parameter optimization searches, pro- 
ducing all cosmological results to within a couple percent of their 'true' values (i.e., the 
values to which the parameters asymptote for much larger pixelizations); and that going to 
~5000 — 10000 pixels yields cosmological output results that are stable to within a small frac- 
tion of a percent. (Though going all the way to ~25000 pixels actually leads to discretization 
problems due to round-off error, which causes problems for higher derivatives of the luminos- 
ity distance function, leading to instability in the value of j® hs ; and so we avoid pixelizations 
this large.) Our high-precision output cosmological parameters - i.e., all numerical results 
other than those from optimization searches - which we will quote in Section [3j below, have 
all been obtained from runs with ~5000 — 10000 pixels (not counting the further addition of 
a smaller number (~1000) of extra pixels used for integrating the Hubble curves out to the 
computationally simpler region past zi nit , before the simulated onset of clumping). 
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3. RESULTS WITH RECURSIVE NONLINEARITIES 

3.1. Weakened Backreaction Effects: Difficulties with "Clumping-Saturated" 

Models 

As mentioned previously, the incorporation of recursive nonlinearities (henceforth RNL) 
does not merely affect the late-time behavior as z — > 0, but in fact exerts a profound damping 
influence upon the entire process of causal backreaction, weakening the overall effect for 
a given set of model input parameters. As an example, we consider the impact of RNL 
upon one of the best-fitting models from BBI, ^Lm(t) with zi nit = 25 and ^ = 1.0. As 
depicted in Figure [2j we see that adjusting the simulations to include RNL greatly reduces 
the apparent acceleration effect - no longer quite even achieving an 'acceleration', in fact, 
since q$ bs = 0.058 > 0. 

For ease of comparison to the results from BBI, in this pap er we will still prima rily 



use the SNIa data from the SCP Unionl supernova compilation (IKowalski et al.l 120081 ) for 



conducting fits of our new models with RNL. As would be expected from Figure [21 the new 
fit to these SNIa data (not shown here), even after re-optimization with respect to H® hs , 
is much worse in this modified model, with its 'chi-squared' value - for 307 SNIa minus 3 
model parameters yielding 304 degrees of freedom - increasing from Xfh = 312.1 without 
RNL, to Xfh = 410.1 with it (where for comparison, xlit = 311.9 for best-fit flat ACDM). 
Visually speaking, though this causal backreaction model with RNL does at least manage 
to clearly separate itself from the strongly decelerating behavior of matter-only flat SCDM, 
it no longer matches best-fit ACDM but now falls significantly short of it, lying somewhere 
in the middle between SCDM and ACDM; and correspondingly, the output cosmological 
parameters - though still being substantially modified from those for the unperturbed case 
without causal backreaction - are now far less reflective of those required for a successful 
concordance, as would indeed be expected from a cosmological model producing too weak 
of an 'acceleration' effect. 

Within the set of 60 model input parameter choices from BBI (all with \l/o < 1), the 
'best' models now, with RNL included, are \I/md with (^i n it,^o) — (5,1-0), which fits the 
SNIa most successfully with x| it = 397.4; and ^sqr, also with (zi nit , \l/o) — (5,1-0), which 
despite yielding a worse overall fit (xfh = 421.3) because of too little clumping at early 
times, manages to produce the largest late-time backreaction effect, achieving an actual 
'acceleration' with g° bs = —0.026 < 0. These two models are plotted in Figure [31 shown 
along with matter-only flat SCDM and ACDM, where all models are now shown in the figure 
as optimized with respect to iJ^ bs for the Unionl SNIa data set. 



One of the most important changes in causal backreaction due to RNL, is the slow- 
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down that it causes for inhomogeneity information propagating in from great cosmological 
distances. This strongly damps the backreaction effects due to perturbations at the outer 
edges of an observer's causal inhomogeneity horizon, which despite being the most distant 
ones are nevertheless very important for backreaction because they include the largest spher- 
ical shells of inhomogeneous matter, containing the greatest incremental mass (per shell) of 
perturbations overall. 

Besides simply weakening the total backreaction effect, RNL acts specifically to inhibit 
the contributions from older perturbations generated during relatively early clumping epochs, 
due to a reduction of the coordinate distance (and thus of the volume of clumped matter now 
enclosed by the sphere) out to some given propagative look-back time; an inhibition that is 
particularly strong if the perturbation potential I(t) has grown too close to unity very early 
on. One important potential result is that this may significantly lessen the ongoing effects 
of causal backreaction into the far future - due in the pre-RNL formalism to eternally- 
expanding observational horizons - an issue related to the ultimate possible fates of the 
universe, to be discussed further in Section I3.4[ below. But our immediate concern here is to 
note that RNL makes later (i.e., more recently generated) clumping far more effective than 
earlier clumping in terms of generating a late-time acceleration. Hence in our runs here with 
RNL, we obtain the strongest apparent acceleration effects with the \I/sqr models, and the 
weakest with the $md models; also, choosing smaller (i.e., later) values of zi nit always makes 
g° bs more 'accelerative' (smaller or more negative), as well. This is all very different from 
the results in BBI, where the \I/sqr functions were by far the weakest models at producing an 
apparent acceleration; and where the $md models were so strong that small or mid-range 
values of zj n i t were necessary to obtain good fits, since those models with the beginning of 
clumping occurring significantly earlier than z\ n \ t ~ 5 led to an actual overkill of acceleration. 

This behavior of causal backreaction with RNL suggests that we may be able to get a 
stronger accelerative effect by using later and later (i.e., smaller) values of 2i nit ; and indeed, 
going from zi n[t = 5 to zj n i t = 2 for the \fsqr model with \& = 1.0 does strengthen the 
observed acceleration from g° bs = —0.026 to g° bs = —0.043. But this type of procedure 
does not save these models as a mechanism for explaining the apparent acceleration, due 
to several obvious problems. First of all, the total amount of generated acceleration is still 
far too small, and thus the fits to the SNIa data remain extremely poor (and in fact, Xmt 
gets worse for the \&sqr runs with later Zj^t)- Besides such problems with the SNIa fits, the 
output cosmological parameters are also unacceptable for such runs, with quantities like t® bs 
and fi™ w growing far too small with decreasing Zi n a, thus spoiling any efforts at achieving 
an alternative cosmic concordance. Lastly, such input parameters do not even make good 
physical sense - can we seriously estimate the effective onset of cosmic clustering as beginning 
at zinit ~ 2, or even later? Clearly, we cannot achieve our goals of reproducing the observed 
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cosmic evolution with causal backreaction in the presence of RNL while stuck within the 
confines of models with ty Q < 1.0. 

These difficulties lead us to the question of what may be wrong with our formalism; a 
question to which there are several possible answers. 

First of all, it is certainly possible that causal backreaction is not the mechanism funda- 
mentally responsible for the observed cosmic acceleration. Even so, given that these effects 
are (as we argue) real and significant - the cosmic expansion is clearly different than it would 
be without them, as we have seen above in Figures [2] and [3], witnessing the clear separa- 
tion between our models and unperturbed SCDM - even a cosmology actually dominated 
by some form of Dark Energy would be strongly modified by causal backreaction, leading 
to misconceptions about the evolution of the universe (and about the nature of that Dark 
Energy, itself) were these effects not taken into account. Alternatively, causal backreaction 
may be part of a 'combination' solution for explaining the apparent acceleration - in con- 
junction, perhaps, with observational effects (recall Section 12. ip such as lensing or other 
inhomogeneity-induced perturbations to the luminosity distance function d^(z) - such that 
the total result from all perturbative effects, combined, may be sufficient to remove the need 
for Dark Energy. 

Another possibility is that our fundamental physical approach towards causal backreac- 
tion is far too oversimplified to properly model the cosmic evolution, just as we have argued 
was true for the standard (non-causal) method traditionally used for estimating backreac- 
tion. There are many possible flaws with our approach - for example, perhaps the linearized 
McVittie solution (from which we calculated our perturbative metric in Section 12.21) is inap- 
propr iate, given that the t otal mass of the perturbation in that metric is fixed as constant in 



time (IKaloper et al.ll2010i ). in contrast to the tendency of overdensities in the real universe 
to keep increasing in mass due to incoming mass flows. (Though our growing \I/(t) functions 
are of course intended as an alternative way of accounting for this.) On the other hand, it 
is quite possible that the complex processes of cosmic structure formation are too dynam- 
ically rich to be modeled accurately by a Newtonianly-perturbed metric containing just a 
single perturbation potential function. The very natures of our \I/(t) functions are obviously 
oversimplified - not accounting at all, for example, for the differences between the simple 
linear evolution of matter overdensities, and the far more complex process of nonlinear den- 
sity perturbation evolution - the latter including subsidiary effects such as virialization and 
stabilization (the detailed time-dependence of which may be quite important, though our 
formalism completely neglects it), as well as "gastrophysics" feedback (e.g., energy injection 
from star formation and from the supernovae, themselves) resulting in baryon shock heating, 
and so on. Considering such complications, it is possible that no simple functional form of 
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any kind for ^(i) - and for that matter, no simple numbers for z^t or \I/o - m &y suffice 
for a realistic description of the physics. Perhaps nothing would be accurate short of a 
fully-detailed cosmic simulation program, calculating a realistically self-consistent represen- 
tation of the intertwined behaviors of cosmic structure formation and astrophysical feedback 
processes as they evolve together over time. 

A different, possibly significant flaw may be our complete neglect of gravitational non- 
linearities; after all, values for the 'Newtonian' perturbation potential in the range of Iq ~ 
0.2 — 0.6 are certainly not negligible, and may imply higher-order gravitational terms from 
general relativity that would contribute significant backreaction effects for any reasonable 
model of the evolution of cosmic structure. Leaving out such nonlinear gravitational terms 
may imply that our calculations here are not serving as an 'estimate' of the overall causal 
backreaction effect, but rather lower limit to it. 

Yet another possibility is that the greatly simplifying assumption of a smoothly-inhomogeneous 
universe is incompatible with a proper treatment of recursive nonlinearities; after all, inho- 
mogeneities massive on a truly cosmological scale, such as galaxy clusters, are not smoothly 
spread in angle, but are concentrated at particular locations on the sky. As such, RNL effects 
would slow down the propagation of inhomogeneity information more strongly along some 
directions than along others; and so perturbative effects from new clumps, located in new 
directions, may be able to 'slip in' with less delay than would inhomogeneity information 
coming in from the same directions on the sky as pre-existing clusters that have already 
imposed their causal backreaction effects upon us. (A kind of cosmic "channeling" effect, 
so to speak.) But quantitatively significant or not, it is likely that such issues cannot be 
properly evaluated without the development of a fully-3D simulation program for causal 
backreaction, as discussed previously in Section 12.11 

On the other hand, the main difficulty here may really be much simpler: the basic phys- 
ical assumptions and methods of our formalism might actually be capable of modeling causal 
backreaction effects well enough, with the problem simply being inadequate choices for our 
(heuristically-adopted) clumping evolution models and model input parameters. Given the 
above results with RNL, it is clear that simply tweaking the detailed time-dependence of the 
\I/(t) functions would not be enough to obtain good data fits; but despite such difficulties, 
there is actually one very simple thing that we can do, which may extend the applicability of 
our causal backreaction simulation formalism, while also dealing with some very real physical 
effects in an empirical yet useful way. From the results in this section, it is clear that we 
cannot achieve a new cosmic concordance with what one might call "dumping-saturated" 
models that are ultimately limited in strength to values of bounded above by unity. But 
what about models with > 1? Such models would clearly lead to much stronger backre- 



-23 - 



action effects (with no pre-determined upper bound, obviously), and could perhaps include 
models that explain all cosmological observations, re-establishing an alternative concordance 
for causal backreaction even in spite of the damping effects of RNL. Given our basic defini- 
tion of as "the fraction of cosmic matter in the clumped, self-stabilized state" , one may 
ask whether it is ever physically appropriate to use values of m excess of unity. In the 
upcoming section, we will argue that the answer is yes, that it is entirely appropriate - and 
in fact, both useful and necessary - for a valid phenomeno logical description of hierarchical 
clustering in the universe. 

3.2. Strengthened Backreaction Effects: Regaining an Observed Acceleration 

with Hierarchical, Nonlinear Clustering 

The convenient but ad-hoc representation of the state of cosmic clumping at some given 
time t with just a single number - ^(t) - is clearly an extreme simplification, as discussed 
above. Even aside from our neglect of the detailed time dependence of virialization processes 
by having this function simply represent the clumped matter fraction over time, there is also 
the fact that large-scale structure is not a featureless blob: the 'clumped' state is not a simple 
binary state, 'on' or 'off'. Most notably, clustering in the universe is not structure- free, but 
instead exhibits a roughly discrete form of hierarchical clustering, with identifiable structures 
existing on a few well-defined scales - for example, galaxies existing as virialized groupings 
of stars, galactic clusters as virialized groupings of galaxies, and so on. This recognizable 
substructure will obviously affect the resulting causal backreaction that is induced, and so the 
question becomes how to best (and most simply) estimate the effects of such substructure. 

One conclusion which we can immediately draw is that clustering with multiple levels 
of substructure will exert more of a total backreaction effect than would a simpler form of 
structure without them. From elementary physical considerations, it is clear that the total 
gravitational potential energy contained within a bound collection of individually-condensed 
clumps will be larger in magnitude (i.e., more negative) than that within a bound volume 
of featureless dust. By the Virial Theorem, the presence of a greater amount of (negative) 
potential energy requires the counterbalance of more kinetic energy, which implies more 
backreaction-causing vorticity and velocity dispersion. Thus the existence of more (stable) 
substructure implies that more induced causal backreaction was generated during the process 
of stabilizing it. 

But one must ask, how much more? Each distinct level of structure actually involves 
the same amount of mass, when considering it over a fixed cosmic volume - whether one 
considers the individual virializations of the tens or hundreds of galaxy clusters in a single 
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supercluster, or of the tens or hundreds of thousands of individual galaxies in that same su- 
percluster, or of the trillions of stellar-mass bodies in that supercluster - and in virializing the 
same amount of mass, it seems reasonable to expect the same amount of volume-generating 
velocity dispersion to occur. Thus we claim that an effective phenomenological way of esti- 
mating the total backreaction effect that has been generated by a stabilized volume of mass 
with multiple levels of internal substructure on different length scales, is by counting the 
effective number of levels of structure that have attained an essentially complete virializa- 
tion. Astrophysically realistic structure formation therefore provides a natural way in which 
ty(t) should be expected to not only reach but to eventually exceed unity, and then to even 
keep growing (without any fundamental theoretical limit) as long as there are new levels of 
structure to be brought into a self-stabilized state. 

To avoid giving ourselves an infinite amount of leeway in searching for an alternative 
concordance, however, it is important to place bounds upon the 'reasonable' values of ^(to) = 
\l/o that one should expect to find when using causal backreaction to explain the apparent 
acceleration. The essence of the question is simple: "How many levels?" In particular, we 
need only consider those typical length scales which became virialized during or not too long 
before the general cosmic time of the observed acceleration (i.e., z < few), since (as will be 
shown in Section [3731) RNL sharply limits the ongoing influence of causal backreaction effects 
from relatively early, 'outdated' virializations (especially when they are strong). 

It seems safe to include at least two scales in our \l/ estimate: that of individual galaxies, 
and that of galaxy clusters, implying a \l/o of at least 2. Going below these classes of structures 
in length scale, there are of course a multitude of individual star clusters and solar systems 
(supported by the vorticity in their orbital motions) that have not too long ago undergone (or 
are still undergoing) relaxation; and even below that in scale, there are continual generations 
of newly-forming solid objects supported by internal body forces (which going beyond the 
perfect-fluid approximation, also contribute positive terms to the Raychaudhuri equation). 
Below this, of course, there is infinite regress; though at some point it becomes absurd to 
consider ever-smaller structures, and significant quantities of new self-stabilization is not 
happening at very small scales in any case. It is difficult to precisely estimate the total 
effective causal backreaction feedback that is induced by structures at sub-galaxy scales; but 
for our purposes, let us consider it a 'significant fraction' of 1. 

At much higher scales, there are groupings of galactic clusters, leading to superclusters, 
and the largest known structures in general. Superclusters are too large to be gravitationally 
bound yet (and certainly were not bound at high redshifts, back around the onset of the 
apparent acceleration), but are experiencing relaxation processes that are underway during 
the current epoch. Once again lacking a precise estimate, we hypothesize that the contribu- 
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tions by stabilized structure formation to backreaction from scales larger than that of typical 
galaxy clusters, in terms of increasing the effective value of may amount to a small but 
non- negligible fraction of 1. 

All told, reasonable values of the present-day clumping evolution parameter for causal 
backreaction, when one considers hierarchical clustering, may be within the range ~ 2 — 3; 
or really pushing it, perhaps within ty ~ 2 — 4. Allowing ourselves such an enhanced 
phase space for fy , the question then becomes whether or not we can re-establish causal 
backreaction, now with RNL included, as a successful (and 'concordant') explanation of the 
observed acceleration. 

For each of the $mdi ^Un, and ^sqr clumping evolution functions, different values of 
Zi n it can be explored, where the highest reasonable value (say, zi n n ~ 25) is set to represent 
a time shortly before the onset of cosmic reionization; and where the smallest reasonable 
value (say, 2i nit ~ 2) would represent the time leading up to serious and intense virialization, 
signaled by such events as the supernovae themselves (noting that the highest-redshift SNIa 
in compilations like those of the SCP Union collaboration tend to be at z ~ 1.5 or so). Then 
for each value of zi nit that we adopt, a search is done over \l/ to find which value achieves a 
best-fit by minimizing Xfh- 

Finding a value of which minimizes Xmt f° r an y & YerL z imt can always be done - and 
often with good results, where the xlit v& l ue is almost as low (and sometimes even lower) than 
that for the best-fit ACDM model. But there is more to the true task of 'optimization' than 
this, since besides achieving a good fit to the SNIa data, one must also match the observed 
values of a variety of complementary cosmological parameters in order to produce a proper 
alternative concordance. Achieving all of these goals simultaneously is more challenging; 
and so to conveniently measure our ability in this regard, we define what can roughly be 
considered an "Average Percent Deviation from Concordance" parameter, as follows: 

1 r |14 RW - 1| |# Obs - 69.96| |# bs - 13.64| 



Avg. % Dev. = _{i_M _L + ^ + 

& 6 1 1 69.96 



13.64 

M»- - (-0.713)1 |/r - 28541 - 311.9) x 

|-0.713| 285.4 311.9 1 ' V ; 

where all of the numerical values used above to represent the 'Concordance' values are 
reference numbers obtained from the Unionl-best-fit flat ACDM model (i.e., Vl\ = 0.713, 
H = 69.96 km s _1 Mpc _1 , with the approximation of zero cosmological radiation density), 
and all quantities are expressed here in typical units (e.g., GYr for t° bs ). 

An optimization search over the range of reasonable input parameters for our three 
cluster evolution functions, using quick simulation runs (i.e., 1000 pixels for t > tinit) , pro- 
duces a number of models that provide both good SNIa fits and good output cosmological 
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parameters. A summary of the results of this search - here showing only those runs which 
possess the optimal ^o.Opt value for minimizing Xfh m eacn particular zi nit case - is given in 
Table [U 

We see that the best (minimized) value of Xmt obtainable is in fact somewhat different for 
different choices of the input parameter z^] but varying even more dramatically with zinit is 
the average deviation from Concordance for its x| it -minimizing run. We have highlighted (in 
bold) five of the cases in Table [TJ ^ M d with (z Init , * ) = (1.5,2.8) and (zinit,*o) = (2,2.9); 
^ Lin with (zyt^o) = (2,3.2) and (z Init , tt„) = (3,3.3); and * Sqr with (z lnit ,V ) = (25,4.1). 
We will refer to these cases as 'best runs', chosen not by the global minimization of Xpit> 
but by selecting zi nit settings for which the x| it -minimizing ^o.Opt run produces cosmological 
parameters with very small deviations from Concordance, in addition to having values of Xpit 
that are acceptably close to the smallest ones found for the relevant \I/(t) function. We will 
consider the astrophysical implications of the particular (z^, ^Opt) input values required 
to obtain these 'best' runs, later on, below; but first, we consider here the principal output 
results of interest. 

These five best runs (now re-done with 10000 pixels each, for higher precision) are all 
good fits to the SNIa data, as can be seen from their Xmt values, which are nearly as small 
as that for the best-fit ACDM model (and are even smaller than that from ACDM for all of 
the \l/sqr runs in the table). The high quality of these fits can also be inferred from Figure HI 
where we plot three of these runs - one each for ^md, ^Um an d ^sqr - against ACDM and 
the matter-only flat SCDM model, showing our runs to be very good at mimicking best-fit 
ACDM (particularly within the redshift region z < 1, containing most of the SNIa), thus 
demonstrating their ability to successfully reproduce the cosmic 'acceleration' as it is actually 
observed. 

The derived output parameters for these five best ^ > 1 runs are presented in Table [2J 
Here we see how their Xfh values, nearly as good or (for ^sqr) better than that for best-fit 
ACDM, makes them comparable to ACDM in terms of the fit probability, Ppit (noting that 
ACDM has fewer optimizable model parameters, and thus more degrees of freedom, giving 
it a higher Pp it value for the same Xmt)- We also see that their present-day 'Newtonian' 
perturbation values, Iq = J RNL (to), do not go much higher than I Q ~ 0.6; this places moderate 
limits upon the effects of gravitational nonlinearities - which are not explicitly dealt with in 
our formalism - implying at least some reasonable level of accuracy for our models regarding 
this approximation. 

Considering the observational cosmological parameters, there is an astrophysically- 
interesting modification of ,2 0bs (relative to z FRW = 1) of ~14 — 17%; while the best-fit 
observable Hubble constant, -£^p bs , is in all cases extremely close (well within 1%) of the 
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ACDM Concordance value of 69.96 km s _1 Mpc _1 . The cosmic ages of these models are all 
at least equal to the "Concordance" value of t° bs = 13.64 GYr (with the ^sqr run being 
exactly equal to it); and since cosmologies with higher (if still reasonable) ages would also 
solve the Cosmic Age Problem just as well, all of these runs (no more than ~0.8 GYr older) 
therefore appear to have quite acceptable values of t® bs . 

The goal of achieving spatial flatness without Dark Energy is also accomplished quite 
successfully, with flf^ ~ 0.94 — 1.16. This result is even more impressive given the fact that 
to produce these numbers, one must first normalize the calculations to some pre-specified 
value of the observed matter density, ^M bs ; an d thus any uncertainties in that observed 
density will translate themselves into error bars on the fiM RW values output from the sim- 
ulations - and with a factor greater than unity, since Af2^ RW ~ Afi^ bs (fi^ RW /f2j^ bs ) ~ 
Afi^ bs (l/f2^ bs ) ~ (3 — 4) x Afi^ bs . For all of the results presented in this paper, we have 
used ^M bs = 0.27. If one assumes that this value is more likely to be too high than too low, 
then models here with results of fi™ W ~ 1-0 would be somewhat better than those with 
^m RW ~ 1-0; but in any case, all of the models in Table [2] are clearly successful at effectively 
achieving flatness. 

In terms of providing a sufficient (apparent) acceleration, these models are again suc- 
cessful, with WQ bs ~ (—0.64) — (—0.75), effectively bracketing the Concordance value of 
w® hs = —0.713. We also note that since t he precise cosmological p arameters obtained from 



SNIa fits are highly model-dependent (e.g.. lCattoen &: Visserll2008l ). it is therefore not neces- 



sary for alternative-cosmology models to exactly reproduce the w^ hs result from Concordance 
ACDM, but merely to generate an 'accelerative' value of u>° bs strong enough to provide a 
sufficiently good fit to the supernova data, as has been done here. 

One final quantity of great importance to be reproduced by an alternative concordance 
is the CMB acoustic scale, Z^ bs - a cosmological parameter that weighs the cosmic evolution 
over vast distances in space and time. We see here that for the Concordance value of 
Z° bs = 285.4 (as computed using a toy ACDM model with no radiation), excellent matches 
are achieved by \&md with (zi n it, 

tf ) = (2,2.9), and ^ Lin with (zi„ it ,#o) = (3,3.3) (which 
is why we chose these $md an d ^Lin models to show in Figure H] over the other highlighted 
ones in Table [fl despite their slightly worse Xfh values); and sufficiently good matches are 
obtained by the remaining three models in Table [2j Although a discrepancy of A/° bs ~ 10 
with respect to the value expected from ACDM Concordance is certainly something that 
would be noticeable from observations, we will see later on (in Section 13. 3j) that one can do 
even better than this with more refined models (for the ^s qr case, particularly). 

The last remaining output parameter quoted in Table |2]is the "jerk" or "jolt" parameter, 
j° bs . As discussed at length in BBI, rather than serving as a well-observed cosmological 
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parameter to be matched by theories aiming for an alternative concordance, the (still poorly- 
measured) value of j° bs can be viewed as a prediction of one's theory, useful for ultimately 
distinguishing such theories from Cosmological Constant Dark Energy (since flat ACDM is 
always constrained to j — 1, for any value of Qa = 1 — fi M ). We recall that the old version of 
our causal backreaction simulation program, without RNL, obtained a set of best-fit models 
which all possessed very high values of the jerk parameter - specifically, jo ~ 2.6 — 3.8 for the 
half-dozen best runs - providing a clear and distinct way of observationally distinguishing 
our causal backreaction paradigm from a Cosmological Constant, or from any other slowly- 
evolving form of Dark Energy reasonably close to ACDM. But here the results in Table [2] 
show that when RNL is properly incorporated into our formalism, the situation becomes far 
less clear: the best-fitting $md and ^Lin models now prefer low values of the jerk parameter 
(jo bs ~ 0.5 — 0.6 and j° bs ~ 0.8 — 0.9, respectively); and the best-fit \I/ Sqr model, though still 
predicting a value higher than unity, now only goes so high as j° bs ~ 1.7. A further analysis 
below in Section 13.31 will demonstrate that even this result may be subject to substantial 
change; and thus the search for a reliable observational method for 'falsifying' our causal 
backreaction paradigm (that is, materially distinguishing it from ACDM) will likely be more 
difficult than what was hoped for in BBI, where the possible impact of recursive nonlinearities 
was still only dealt with in terms of qualitative caveats. 

With this analysis of our 'output' results done, we now look at the input parameters 
which were necessary to produce these good cosmological fits. First of all, considering z Init , 
we see from Tabled] that to obtain good cosmological output parameters from the $md an d 
^Lin models, we have had to impose an extremely late value for the time representing the 
'beginning' of structure formation - specifically, z^a ~ 1.5 — 3. This is a direct result of 
the effects of RNL, which for an already-perturbed universe will slow down the continued 
propagation of old perturbations coming in from great cosmological distances, thus rendering 
the causal backreaction contributions of early-developing inhomogeneities to the late-time 
'acceleration' much weaker than they were for the equivalent models simulated in BBI. This 
is a serious problem for the and (even more so) for the $md functions, which are 
especially designed to provide more early clustering, at the expense of late-time clustering. 
We must ask, in fact, whether input values like these for z Init even make astrophysical sense, 
given the fact that z° hs ~ 3 actually marks the epoch during which the large-scale collapse 
of material and the resulting galactic feedback (from star formation, etc.) begins acting 
strongly to shock-heat c osmic baryons, sendin g a significant portion of the baryons into the 



superheated IGM (e.g., ICen fc Ostrikerll2006l ). thus helping to slow down further structure 
formation. 



This problem goes away, however, for the \l/sqr runs, which have enough late-time clump- 
ing to remove the need for choosing a late (i.e., small) value of Zi n a to get good output cos- 
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mo logical parameters. The runs only get better, in fact, as one goes to higher Zi n a, and we 
stop at Zi n it = 25 not because optimization tells us to (though increasingly early initialization 
times eventually brings diminishing returns), but because it is not physically reasonable to 
set the start of clustering too long before the onset of reionization. 

Besides this success of providing an alternative concordance with reasonable zi n n values, 
some reflection shows that the ^sqv runs are the most sensible in every regard. As we now 
consider the *ff(t) function to represent the feedback effects of completed clustering on several 
different, dynamically-defined scales in hierarchical structure formation - rather than as a 
simple percentage of 'clumped' versus 'unclumped' matter - the causal backreaction process 
becomes more naturally focused upon structure formation in the highly nonlinear regime of 
density perturbations. And as we recall from above, density fluctuations in the nonlinear 
regime evolve as 5p/p oc t m with m > 2 for a matter-dominated universe, meaning that 
ty{t) functions evolving as \l/sqr oc t 2 (or functions growing with even larger powers of t) are 
therefore the most relevant models for causal backreaction, rather than those with smaller 
powers of t like the ^T.in or $md models (the latter type evolving as Sp/ p in the linear regime), 
which had been the dominant ones for generating backreaction in the old simulations without 
RNL. 

Considering the above results once more in this light, we see that \l/sqr has indeed 
performed better than or ^md, for example looking much more like best-fit ACDM 
in Figure HI thus providing better SNIa fits with significantly lower Xmt vam es (recall Ta- 
bles [Tj2])- Furthermore, the output cosmological parameters from the \l/cropti m ized \[^ Sqr runs 
(with sufficiently large zi nit values) were all extremely good, with just the possible exceptions 
of a slightly low total density (e.g., f^ RW ~ 0.94 for zi n \ t = 25), and a slightly errant CMB 
ane ular scale (A/£ bs ~ 9); but even those discrepancies are quite small given the very simpli- 
fied nature of these clumping evolution models. (And those discrepancies can be made even 
smaller with some further development of the models, as will be seen shortly, in Section l3TB~j) . 

The only real problem with these \l/sqr models, and quite an obvious one, is the very 
high degree of clumping needed in order to generate these results - that is, ^o.Opt — 4.1. 
This is right around the extreme upper edge of the 'reasonable' range of values for ty , 
which we previously argued should likely lie within ty ~ 2 — 4 (and more conservatively, 
within $0 ~ 2 — 3). Such \l/o,Opt values for \l/s qr (which are greater than those for or 
$md) ar e large enough as to potentially strain the credibility of causal backreaction as a 
cosmological solution. As we will see below, however, this concern can be alleviated simply 
by adding in one important and physically reasonable model parameter as part of these 
clumping evolution functions: specifically, a parameter representing information about the 
possible end of clumping, as could be triggered in the real universe by galactic feedback. 
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3.3. Early-Saturation Clustering Models: Parameter Flexibility and Paradigm 

Falsifiability 

As our clumping evolution models have so far been denned, each takes on its 'ultimate' 
value, ^o, at t ; and this may make physical sense in a universe where the effective dumpiness 
is still growing, so that the value ^/(t ) is not particularly special, but just represents a 
'snapshot' of ^(t) taken at our particular observation time. But if there is instead some halt 
to the growth of clumping - due either to the limiting value of \I/(t) — » 1 for our pre-RNL 
simulations, or due to various feedback effects in these new models (which no longer possess 
a pre-determined upper bound on ^(t)) - then from a Copernican view, there is nothing 
necessarily special about the current cosmic epoch, and so there is no reason to assume that 
the asymptotic value of ^(t) (if there is one) should occur specifically at the present time. 

Having already re-evaluated our formalism by interpreting the implementation of \l/o > 1 
as representing multiple levels of structure on different cosmic scales, it has become possible 
for ^(t) to continue increasing without any hard theoretical limit for t > t , perhaps con- 
siderably beyond the effective degree of clumping that exists today. On the other hand, if 
there does happen to be some naturally-occurring, effective limit of clumping, then there is 
no reason to believe that the processes in charge of cosmic clustering would wait until t — > to 
to reach it. 



As noted above, ICen fc Ostrikerl (120061 ) show that the nonlinear collapse of material and 
feedback from star formation have likely acted to shock-heat cosmic baryons to millions of 
degrees, thus inhibiting clumping - the baryonic part of it, at the very least - by keeping 
and/or sending a significant portion of the baryons into the superheated IGM; their Figure 
1 shows simulation results depicting the rapid increase of very hot baryons in the IGM 
beginning strongly as of z° hs ~ 3. 

Additionally, the backreaction effect itself may be helping to put the brakes on structure 
formation at late times; particularly given the fact that it is generated most strongly within 
the vicinity of the virializing masses themselves, in high-density regio ns. In that case, signs 



of a cosmologically recent slow-down in the growth of clustering (e.g.. IVikhlinin et al.ll2009[ ) 
- findings which are viewed as strong evidence in favor of Dark Energy - may instead 
be indicating the effects of causal backreaction. (Conceivably, feedback from backreaction 
may also be a contributor to other poorly understood cosmological effects such as galaxy 
downsizing ( Cowie et al. 19961), the cuspy CDM halo problem, and the possible dearth of 



dwarf satellite galaxies (e.g.. iPrimackl 120031 ) . 



What all of this implies is that our clumping evolution functions, as they have been 
defined in Equations [T^h-c. would likely not remain meaningful all the way to t — > to; and 
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thus it is useful to consider how the functional form of our models, ^(i), might be altered to 
account in some manner for this effect, to represent the possible late-time softening of the 
growth of clustering due to astrophysical and/or backreaction-related feedback processes. 

The most straightforward way to model this behavior is by implementing an early "sat- 
uration time", tsat < (with ^Sat = [(^oAsat) 2 ^ 3 — 1]), representing the time before the 
present epoch at which point the growth of clustering became saturated to a final, fixed 
value, ^(tsat) = ^o- This new clumping evolution function then fixes \l/(t) to remain equal 
to \&o ah the way from t Sat to t . Recalling from Section I3T21 that the $md and ^un models 
only obtain their best fits - and their most 'concordant' cosmo logical parameters - for very 
late ti n it values (i.e., zi n i t ~ 1.5 — 3), it then makes little sense to outfit those models with 
an even later tsat > Wtj we therefore make this modification only for the \l/sqr model, which 
was optimized for concordance with much earlier t Init values. The resulting, retrofitted \&sqr 
function for early saturation is thus defined as: 



This is obviously a very elementary way of accounting for how astrophysical feedback acts to 
slow down the growth of clustering (with the resultant lessening of its causal backreaction); 
but even such an oversimplified model will provide us with valuable cosmological lessons. 

As a first run to explore the implications of early saturation, we consider again our 
earlier, cosmologically successful run with \&sqr (functionally identical to \I/sat with ^sat = 0), 
which had used zi nit = 25 and its corresponding \I/o,Opt = 4.1; these results are now compared 
to those from a new \l/sat run still using 2 Init = 25 and (the now non-optimal) = 4.1, but 
with z SsLt = 3 now imposed - a conservatively early value suggested by the onset of substantial 
hot baryon injection into the IGM. 

The resulting residual Hubble diagrams for the two cases are compared in Figure 
where the powerful impact of moving from zs a t = to ^sat = 3 upon the time-dependence of 
the observed backreaction can clearly be seen. Not unexpected is the strongly-amplified ap- 
parent acceleration effect at early times (z° hs ~ 10), due to the greatly compressed timescale 
over which clustering grows from ^(t) = to *&(t) = \I/o = 4.1. But what may be surprising 
is how dramatically the effects of these early-developing inhomogeneities become virtually 
irrelevant at late times, with the effective 'acceleration' almost completely fading out by 
2; 0bs ~ 1 — 2, and the residual Hubble diagram relapsing to a nearly perfect SCDM (i.e., 
decelerating) cosmology not long after z ~ 1. 

This is in sharp contrast to what would be expected from the results explored in BBI, 




(21) 



-32 - 



where early-developing clumping provided the strongest ongoing causal backreaction effects 
due to the ever-growing "inhomogeneity horizons" of observers. The fact that this is no 
longer the case is clearly due to RNL, which damps such ongoing backreaction effects by 
slowing down the causal propagation of inhomogeneity information (as well as by diluting 
the perturbation effects of already-seen inhomogeneities via the extra volumetric expansion), 
thus greatly restricting the cosmological horizon out from which developing perturbations can 
affect the observer. As can be inferred from Equations [T0][T7l above, dr/dt oc a/1 — J RNL (t) 
for inhomogeneity information propagating towards an observer at null-ray speed; and so 
models with very strong early clumping - particularly those with a very early Zsat ~ drive 
up the value of J RNL (t) so close to unity, so early on, that the propagation of inhomogeneity 
information is practically frozen to a halt at later times. The cessation of new inhomogeneity 
information reaching the observer leads to a corresponding freeze in the continued evolution 
of 7 RNL (t); and as is obvious from our metric given by Equation [TJ a constant J RNL value 
can simply be transformed away via redefinitions of t and r, thus making a cosmology 
with static I(t) = J RNL look exactly like a decelerating, matter-dominated SCDM universe. 
(Interestingly, a considerably stronger late-time acceleration effect can be generated by using 
a smaller value of \l/o for these very early zsat runs, which lessens this 'freezing' effect. The 
best possible SNIa fit for this zsat = 3 case is therefore achieved with the relatively low value 
of \I/o,Opt = 1-2; though even that remains a very poor fit in absolute terms.) 

The importance of such behavior is that it limits the degree to which causal backreaction 
with RNL can be 'self-powered' via the ever-expanding reach of observational horizons, in 
defiance of a static final value; a result which increases its dependence upon being 'driven' 
by a continually-growing ^(t) function. As Figure [5] has shown, an early saturation of 
clumping can lead to a virtually complete shutdown of apparent acceleration not much later. 
(Though noting again the caveat that nonlinear gravitational effects are not accounted for 
in any of our models - and that they may indeed be important here, since all of the very- 
early saturation models with zsat = 3 which we tested do have excessively large Iq values, 
all significantly exceeding 0.7, with some nearly approaching unity.) These considerations 
will have important implications for the eventual fate of the universe, as will be discussed 
below in Section I3.4( but first, we consider their ramifications regarding the past cosmic 
'acceleration' that has already been observed via the supernova data sets. 

Since a value of zsax = 3 effectively terminates the impact of causal backreaction far 
too prematurely for it to account for the apparent acceleration seen in the SNIa data from 
around z° hs ~ 1 to now, we consider runs with ^sat values that are similar to or smaller 
(i.e., later) than this epoch of acceleration. Specifically, we choose values of z Sa _ t = 2 Iat W = 
(1.0, 0.5, 0.25) for study - i.e., t|at W — (0.35, 0.5, 0.7) - representing times by which the mass 
fraction of cosmic baryons possessing temperatures of T > 10 5 K has increased to reach 
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0.4 — 0.5 (jCen fc Ostrikerl 120061. Figure lb). These values of zsat are also well-placed to 



further illuminate the results of IVikhlinin et al.l ( 120091 ). which demonstrated the significant 



alteration in clustering behavior (due either to Dark Energy or to causal backreaction) 
between one cluster sample with z > 0.35 (and particularly its most distant subsample with 
z > 0.55), versus their more nearby cluster sample at z < 0.25. 

Re-optimizing \l/ for each of these chosen ^sat values (with Z[ nit = 25 fixed for all of 
these runs), the resulting 'best runs' with early saturation are those with the parameters: 
(zsa*,*o) = (0.25,2.6), («sat,*o) = (0.5,2.3), and (s Sat ,¥ ) = (1-0,2.2); we compare this to 
the best-fitting zs at = run, which had ^o.Opt =4.1. 

Residual Hubble curves for these four runs are plotted in Figure O where we see that 
while these early-saturation runs are no longer effectively indistinguishable from best-fit 
ACDM (as was the zg a t = run), they are nevertheless very close to it - particularly within 
the redshift range containing most of the SNIa - thus providing quite good fits to these data. 

The complete fit quality results and cosmological output parameters for these runs are 
presented in Table El First, we note that there are limits to the largest value of Zsa,t which 
may practically be used, since we see that as zsat gets as high as 1.0, the fit probability 
Ppit gets progressively worse (i.e., xfit begins to grow unacceptably large), the cosmological 
parameters become increasingly 'discordant' (particularly fi^ RW and /A bs )' an d there are 
even concerns that the formalism itself begins to break down due to gravitationally nonlinear 
effects, given the disturbingly large value of Iq. 

If we restrict ourselves to ,2sat < 0.5, however, then at a cost of only a small increase 
in Xfhj the situation gets considerably better with the use of nonzero Zs&t in several crucial 
ways. Considering the ^sat = 0.25 case in particular, moving away from z Sat = has improved 
the match to the CMB data and the verification of spatial flatness (better Z° bs and f2jy[ RW 
values, respectively), while remaining essentially as good in terms of the other cosmological 
parameters, with just a tiny degradation in the fit probability (by APpit = —0.025). But 
most importantly, we see that all of this been achieved with a much lower \l/o value - Table [3] 
showing how \l/o,Opt decreases with increasing zs at , in general - dropping all the way from 
^o,Opt = 4.1 for z Sat = 0, to ^o.Opt = 2.6 for z$ at = 0.25. This is now well within the 
\I/ ~ 2 — 3 range of 'reasonable' values for hierarchical clustering, as was specified earlier in 
Section [321 

(In addition, we have also conducted fits to a m ore recent supernova data set, the 



SCP Union2 SNIa compilation ( lAmanullah et al.l 120101 ); and when we minimize Xmt f° r the 



zsat = 0.25 case with respect to the Union2 data, rather than the Unionl data, we still get 
similar results: an excellent fit with Xmt = 

543.6 (compared to xht = 542.7 for ACDM) 
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and Ppit — 0.6; only a slightly higher clumping strength value of \l/o,Opt = 2.8; and similar 
cosmological parameters that are also quite acceptable. Only the unperturbed matter density 
is slightly high with Qf^ w = 1.13; but as discussed above, this may be due to the adopted 
value of fi^ bs = 0.27 being slightly too large.) 

In consequence, it is justified to say that we have in every sense re-obtained a successful 
alternative concordance with causal backreaction in the presence of recursive nonlinearities, 
having generated the proper amount (and temporal behavior) of apparent acceleration - 
as evidenced by a fit to the SNIa data that is very nearly as good as that achievable with 
best-fit ACDM (even without performing a rigorous x Fit -minimization search over our model 
parameter space) - as well as having produced output cosmological parameters that are more 
than acceptably consistent with a variety of complementary cosmic measurements. The price 
to be paid to achieve this goal, due to the incorporation of recursive nonlinearities, is the 
necessity of permitting \l/o values in excess of unity - which though going against the original 
interpretation of this (heuristic) model parameter, as defined in BBI, does seem justifiable in 
a realistic cosmology exhibiting hierarchical structure formation on a variety of scales. On 
the other hand, a new benefit of these results is the shifted focus from f MD models to ^sqr 
models, equivalent to a shift in emphasis from linearized matter density perturbations to 
nonlinear density perturbations; a new emphasis which in fact makes much more sense for a 
causal backreaction paradigm depending upon vorticity- and velocity-dispersion-generated 
virialization as the fundamental origin of structure-induced perturbations to the observed 
cosmological metric. 

One last (yet very important) consideration still remains regarding early saturation, 
though, relating to its impact upon the very late-time behavior of the cosmic evolution. One 
way to approach this, is to consider why the optimized value of ^ drops so precipitously with 
increasing z Sat , so conveniently solving our problem by reducing ^o.Opt to believable values. 
To understand why this happens, consider the clumping evolution functions, themselves, for 
the early saturation models that we have studied; plots of these \&(t) functions are shown in 
Figure 

From these plotted curves, we see that despite the vast differences in ^o,Opt f° r the 
z Sa,t = case versus those with z Sat ^ 0, as long as one does not use too large a value of 
2sat (sticking to, say, ^sat < 0.5), then the degree of clumping ty(t) at mid-range values of z 
(e.g., z FRW ~ 0.2 — 0.5) actually remains fairly similar from run to run. The distinctively 
huge increase from ty(t) ~ 2.5 to \I/(t) ~ 4 for the zs&t = case does not actually happen 
until very late times, ^ FRW < 0.2. But why should completely different behaviors of \l/(t) at 
late times have such a small effect upon the fits to the SNIa data? The answer is twofold: 
first, due to a lack of sensitivity of the data to such changes; and second, due to the nature 
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of causal backreaction, itself. 

For the first issue, consider that most of the SNIa included in these supernova compila- 
tions for measuring the cosmic acceleration are located at fairly high redshift; it is therefore 
unsurprising that the evolution of ^{t) after the epoch of these supernovae should have lit- 
tle effect upon the acceleration measured by those high- z SNIa, regardless of whether ^(t) 
experiences a continued increase or a saturation. But, even if lower- z SNIa were also in- 
cluded in the study, and could be trusted to accurately map out the Hubble flow despite 
their local peculiar motions, there should still be little effect other than a late-time offset, 
perhaps registered as a small change in if^ bs . As noted by iLinderl (120071 ). simple cosmic 



evolution functions that are relatively insensitive to time variations tend to measure an av- 
eraged cosmological equation of state around a 'pivot' redshift of about z° hs w 0.4. And so 
it is understandable why the detailed behavior of after this crucial epoch should end 
up having little effect upon the precise amount of apparent acceleration measured. 

For the second issue, we recall from Equation [5] that the amount of causal backreaction 
due to a spherical shell at coordinate radius r will be proportional to dM/r oc r 2 jr oc r, and 
so very late clumping - which corresponds to relatively small look-back times t ret , and thus 
small distances r - will simply not involve a sufficiently large amount of inhomogeneous mass 
to generate significant causal backreaction. Thus the detailed dynamics of inhomogeneities 
at higher z (and larger r) - assuming distances not so far away as to get damped by RNL, or 
having look-back times so far back relative to ti nit as to have very small clustering, ^(t rct ) - 
will be more important in determining the detailed effects of causal backreaction than would 
the very late (z — > 0) behavior of clumping. 

The overall result is that there is a parameter degeneracy for causal backreaction models, 
where one can reduce \l/o i n tandem with an increase in zgat> without much change in the 
quality of the fit to the actual SNIa data. This degeneracy has a beneficial aspect, in that it 
has allowed us to successfully create an alternative concordance with causal backreaction - 
and to do so in a flat, matter-only universe without any form of Dark Energy - while using 
astrophysically realistic values of \I/o- But it has a negative aspect, as well, in that it lessens 
our ability to predict precise ranges for the late-time cosmological parameters that would 
be output by the various (successfully concordant) models in this formalism. Measurable 
cosmological parameters such as i?° bs , w® hs , and jQ hs are defined - in theory, at least - in 
terms of the z° hs — > behavior of Taylor expansions of the luminosity distance function, 
d^(z 0hs ). Hence, while a variety of ^/(t) functions may all provide good fits to the key SNIa 
located at mid-to-high redshifts, these functions may all have quite different t — > behaviors, 
and thus very significant differences in their cosmological output parameters. 



In particular, consider again the first three of the four runs in Table [3j As one goes 
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from zsat = to zs a t = 0.5, we see that H° hs decreases below the ACDM Concordance 
value of HQ hs ~ 70; but yet, the difference is small enough to not be a serious concern. 
For WQ hs , obtained from one more differentiation of g?l, the difference is larger, going from 
—0.75 to WQ hs ~ —0.59; and yet, the precise numerical value of the cosmic 'equation 
of state' WQ hs is of no great concern to us here - since it is not our task to pin down the 
physics of some hypothesized form of Dark Energy, but merely to reproduce the cosmological 
observations - so long as a good fit with a sufficient amount of apparent acceleration can be 
produced, as has indeed been accomplished. 

The biggest change, however, occurs for j° bs - obtained via yet another differentiation 
of c?l - which drops all the way from j° bs ~ 1.7 for zsat = 0, down through the ACDM value 
of j = 1 and far past it, even going negative for z Sat values that still provide decent SNIa 
fits. 

A small part of this change may be due to the highly simplified nature of the ^sat 
functions that we use; in particular, since we obtain a Hubble curve through two integrations 
of ty(t) - that is, -»■ I(t), and then I(t) ->■ d L (cf. Equations loT TTTj) - the j° bs parameter 
(obtained from three differentiations of d^) essentially contains one differentiation of ^{t). 
As is obvious from Figured our simple ^sat functions are not differentiable at zsat, meaning 
that q° hs (z) (and thus w (z)) will also be non-differentiable there, and so j° hs (z) will pick 
up an actual discontinuity there. Still, the discontinuity is not large enough to account for 
the huge overall differences in j^ bs for the different Zsat values; and most of the effect of the 
jump in j 0bs (z) at Zgat is almost certainly real, due to the real change in the cosmic evolution 
at that point, and would likely still take place to a very similar degree (just more spread out 
in time) for some more complicated \I/sat functions designed to apply smoothing at zgat- In 
any case, the set of runs in Table [3] does fairly clearly establish the trend that increasing 2; Sat 
away from zero leads to a steady decrease in the final value of j^ bs . 

What we can conclude with some certainty, unfortunately, is that we have lost any rea- 
sonable degree of predictability for this j° bs parameter in the causal backreaction paradigm, 
due the weak dependence of the backreaction upon the very-late-time behavior of ^f(t), and 
the resulting degeneracy in (zsat^o)- This is significant, because of the different possible 
ways suggested in BBI to distinguish between causal backreaction and Cosmological Con- 
stant ACDM (or any similar form of Dark Energy), in order to provide our paradigm with a 
falsifiable test, the clearest signature by far was the search for j^ bs ^> 1. For all intents and 
purposes, the reliability of this signature now appears to be gone, and some more intricate 
means will obviously be needed to distinguish causal backreaction from even the simplest, 
A-like version of Dark Energy (though j° bs ^ 1 would still rule out ACDM in favor of some 
alternative model). Thus while the use of early saturation has greatly enhanced the prospects 
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for achieving an alternative cosmic concordance with realistic values of the clumping param- 
eter \l/ ) the incorporation of this additional degree of freedom in the models has made it 
significantly more challenging to definitively replace the paradigm of Dark Energy with that 
of causal backreaction, without resorting to aesthetic arguments of a subjective nature. 

3.4. The New Cosmic Future: How Recursive Nonlinearities Affect the 
Possibility of Long-Term Acceleration from Causal Backreaction 

If appropriate observational tests should eventually be able to demonstrate causal back- 
reaction as superior to Dark Energy as the driver of the cosmic evolution, then one final 
question of great importance would of course be: What is the ultimate fate of the universe? 

This topic was discussed in detail in BBI, in which it was found that the long-term 
cosmic fate depends upon which way the balance tips between the forces working to power 
the cosmic 'acceleration', versus the opposing factors acting to restrain it. For the former, the 
only real influence working to keep the apparent acceleration going (and perhaps ultimately 
promote it in strength to a real volumetric acceleration) was the ever-expanding causal 
horizon, growing in time, out from which an observer can 'see' substantial inhomogeneities 
- i.e., the farthest distance from the observer out to which t iet > tj nit , ^(t ret ) 3> are still 
true. 

For opposition to the acceleration, one source of restraint is the FRW expansion of the 
universe itself, which exerts a natural damping effect upon causal backreaction by diluting 
the density of inhomogeneities, simply by pulling them farther away from one another (and 
from any given cosmological observer) over time. Also acting to limit a possible long-term 
acceleration was the inevitable saturation presumed for ^(t), which - in our original inter- 
pretation of this measure of clumping, as described in BBI - would be limited by an upper 
bound of * ->■ 1. 

Considering all of these factors, analytical approximations were derived for the (pre- 
RNL) metric perturbation function I(t) as t — > oo; and it was found that I(t) always 
asymptotes to a constant numerical value for \I/(t) functions evolving as a simple power of 
t, ^(t) oc t N . This asymptotic value of I(t) is larger for smaller N, with it being equal to 
unity (implying a complete breakdown of the Newtonianly-perturbed metric) for N = 2/3 
(i.e., ^ = $md)- Thus a fully general-relativist ic acceleration at late times due to causal 
backreaction did appear to be realistically possible, according to the formalism of BBI - 
though of course that conclusion was made pending the still-undetermined effects of recursive 
nonlinearities, and other complications. 
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In this paper, however, we will have to revise our expectations, since the incorporation 
of RNL (recursive nonlinearities) has made the prospects of an 'eternal', self-powered ac- 
celeration seem far less likely. We can conclude this despite the fact that I(t) is no longer 
necessarily amenable to a simple analytical analysis, even for t — > oo. One point in favor of 
such a conclusion is that RNL has forced us to change to ^(t) oc t N functions with larger 
exponents N in order to re-establish a concordance - i.e., moving from $md and ^tia before, 
to ^sqr now ~~ an d as our previous results have shown, larger exponents in *ff(t) lead to 
smaller late-time values of I(t) (even if we no longer know its exact asymptotic behavior), 
thus implying a less general-relativistic, less 'accelerative' long-term evolution. But an even 
more important factor is what we learned about the impact of the ultimate 'saturation' of 
the clumping evolution function at some final, fixed numerical value, as was depicted 
in Figure This result (and other simulation runs that we have done along these lines) 
clearly demonstrate that RNL, which acts to stall the expansion of an observer's "inhomo- 
geneity horizon" - particularly so in cosmologies with strong early-time backreaction - will 
largely 'shut off all apparent acceleration effects not long after the ongoing clustering and 
virialization have ceased, once has settled down to some mostly- const ant final value. 
This behavior certainly works against the possibility of an 'eternally' (or even long-term) 
'accelerating' universe. 

The one enhancement introduced in this paper which might possibly help lead to a 
long-term acceleration, is the now fundamentally unbounded nature of ^(t) due to hierar- 
chical clustering (assuming that an effective end to clustering is not imposed explicitly by 
early saturation, as we did in fact choose to impose for our models in Section 13.31) . But 
while this change in interpretation allowing > 1 is indeed based upon real physics, 
is it strong enough to keep the 'acceleration' going continually, deep into the future? The 
answer to this question is naturally uncertain: on (relatively) small scales, clustering has 
never ceased, as star clusters and small galaxies continue to merge (with new virialization) 
into large galaxies; galaxies continue to merge into galactic clusters; and so on. Yet, one 
would realistically expect to find steadily diminishing returns on such 'small' scales. Any 
real hope for a continuing acceleration would seem to rest upon the future clustering be- 
havior on extraordinarily large scales - a realm of structure formation with no real upper 
limit, as superclusters eventually manage to self-virialize internally, then begin themselves 
to merge into even larger structures, and so on, ad infinitum. Higher and higher scales of 
clustering take exponentially longer to complete than the levels below them, though, and 
it is not clear how large the effective 'inhomogeneity horizon' for causal backreaction can 
ever practically become - especially given the fact that the cosmic acceleration itself tends 
to impede the expansion of observational horizons. A truly long-term cosmic acceleration 
(apparent or real) into the far future would therefore seem very difficult to accomplish using 
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causal backreaction with recursive nonlinearities; though such conclusions cannot be con- 
sidered definitive without more detailed cosmological simulations, employing a significantly 
more sophisticated treatment of the effects of vorticity and virialization than we have used 
in these calculations so far, given our toy models for *&(t), and therefore I(t). 

There are still two remaining wildcards, however, which were discussed in BBI and must 
be mentioned again here. The first one is the possible advent of truly general-relativistic 
perturbations, causing the breakdown of our Newtonianly-perturbed metric approximation 
(Equation [2]), and the failure of our treatment of the individual perturbations due to the 
multitude of locally- clumped masses as being linearly summable. The possibility of such a 
complication still exists, as is made obvious (for example) by the very large Iq values that we 
typically find for early-saturation runs with large zg at settings (e.g., the ^sat = 1 case from 
Table [3]). As noted above in Section [3.31 the fact that the propagation of inhomogeneity 
information obeys the proportionality dr/dt oc yL — J RNL (t), means that the flow of such 
information will be choked off whenever J RNL approaches unity, thus freezing the expansion 
of the relevant inhomogeneity horizons, thereby locking I RNL (t) nearly in place at whatever 
actual value that it has managed to grow to by that time. (This is in fact the biggest 
reason why the zsat = 3 model from Figure [5] experienced such an abrupt shut-off of its 
apparent acceleration so quickly after the growth of its clumping evolution function \&(t) 
had ceased.) What we cannot know from our formalism, however, is how to physically 
interpret the effects of a metric perturbation function J RNL (/f:) that is hovering around a 
nearly fixed value, but where that value is always slowly asymptoting towards unity. Does it 
look almost exactly like a decelerating SCDM universe, since a constant I RNL (t) can simply 
be transformed away through coordinate redefinitions? Or would it instead be a universe 
perpetually appearing to be right on the verge of undergoing a runaway acceleration, just 
never quite doing it, basically 'riding the edge' forever? (Perhaps succeeding, at least, at 
producing a hesitating, stop-and-go accelerating behavior.) Or on the other hand, does 
the proximity of 7 RNL (t) to unity manage to override all other considerations, and lead to 
an actual, volumetric, runaway acceleration? The real answers to these questions cannot 
be determined without at least some nonlinear gravitational terms being added into the 
formalism, if not actually requiring a fully general-relativistic treatment. But while the 
quantitative calculation of such effects to determine their physical implications is beyond 
the scope of the causal backreaction formalism that has been presented here (or in BBI), it 
is clear that the sum of innumerably many Newtonian-level perturbations is indeed capable 
of driving the total cosmological metric perturbation right up to the breaking point, where 
a 'real' cosmic acceleration (by any definition) may very possibly take over, and for an 
indeterminate period of time. 

Last of all considerations, though, is the one that can least safely be neglected: and 
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that is the eventual breakdown of our smoothly-inhomogeneous approximation, as increasing 
dumpiness leads to fundamentally non-ne gligible anisotrop i cs acr oss much of the observable 
universe itself (e.g., the "Dark Flow" of iKashlinsky et al.l 120101 ); a scenario which in BBI 
was termed the "Big Mess". This final, virtually certain breakdown of the Cosmological 
Principle is the ultimate game-changer, and questions about the possibility of a 'permanent' 
cosmic acceleration due to causal backreaction then become just as hard to define as they 
are to answer, as the fundamental FRW basis of cosmological analysis finally breaks down 
entirely. 



4. SUMMARY AND CONCLUSIONS 



In this paper, we have revisited the causal backreaction paradigm introduced in lBochner 



( 120111 ). for which the apparent cosmic acceleration is generated not by any form of Dark En- 
ergy, but by the causal flow of information coming in towards a typical cosmological observer 
from a multitude of Newtonian-strength perturbations, each one due to a locally clumped, 
virializing system. Self-stabilized by vorticity and/or velocity dispersion, such perturbations 
are capable of generating positive volume expansion despite their individually-Newtonian 
natures. Noting that previous 'no-go' arguments against Newtonian-level backreaction are 
based upon non-causal backreaction frameworks, we see that the sum total of these small 
but innumerable perturbations adds up to an overall effect that is strong enough to explain 
the apparent acceleration as detected by Type la supernovae, as well as permitting the 
formulation of an alternative cosmic concordance for a matter-only, spatially-flat universe. 

Our purpose here has been to develop and test a second-generation version of this causal 
backreaction formalism, filling in one of the most important gaps of the original 'toy model' 
by including what we have termed "recursive nonlinearities" - specifically referring to the 
process by which old metric perturbation information tends to slow down the causal propa- 
gation of all future inhomogeneity information, therefore reducing the effective cosmological 
range of causal backreaction effects, and thus damping the strength of their overall impact 
upon the cosmic evolution and upon important cosmological observations. 

Utilizing the new simulation program introduced here, which now incorporates recursive 
nonlinearities into causal backreaction, we find profound differences in the resulting cosmo- 
logical model calculations. For a given magnitude of self-stabilized clustering assumed for 
large-scale structure, denoted by dimensionless model input parameter \l>o, the overall power 
of causal backreaction is now considerably weaker, in addition to fading out relatively rapidly 
after the growth of clustering ceases. This is unlike the results of the original model, in which 
causal backreaction effects would continue to grow regardless of any late-time saturation of 
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clustering, due to the causally-expanding "inhomogeneity horizon" seen by an observer which 
continually brings more 'old' inhomogeneities into view from ever-greater cosmic distances. 

After discussion of some of the possible reasons for which causal backreaction may now 
appear to fall short of its cosmological goals - either due to issues regarding the fundamental 
mechanism itself, or due to our highly simplified treatment of it - we then considered a 
very straightforward way in which the paradigm may be fully recovered as a cosmological 
replacement for Dark Energy: all that is needed is the adoption of \l/o values greater than 
unity. Though representing an ad-hoc modification of the original formalism, the change 
makes astrophysical sense in a number of ways. Rather than viewing the clumping evolution 
function ^(t) as simply representing the fraction of cosmic matter in the 'clumped' versus 
'unclumped' state at any given time, ^(t) can now be recognized (more realistically) as 
representing the total backreaction effect of hierarchical structure formation in the universe, 
where clustering and virialization take place simultaneously on a number of different cosmic 
length scales - from stellar clusters, to individual galaxies, to galaxy clusters, etc. Model 
input parameter \l/o = ^(^o) is now interpreted as the effective number of 'levels' of completed 
clustering that exists (at current time to) i n the large-scale structure when one sums over 
the (partial or total) clustering of matter on all relevant cosmic scales. 

Given this enlarged parameter space with \I/o > 1 now permitted, we once again find a 
selection of model cosmologies that succeed (despite the damping effects of recursive non- 
linearities) at reproducing the observed cosmic acceleration, while also re-establishing an 
alternative cosmic concordance by producing output parameters that match the observables 
derived from several of the most important cosmological data sets. Furthermore, astrophys- 
ical considerations regarding the necessary input parameters for these apparently successful 
models - specifically, the need to assume a sufficiently early beginning of clustering - results 
in a preference by the new formalism for ty(t) models that reflect the late, nonlinear phase 
of structure formation. This is an improvement over the old formalism without recursive 
nonlinearities - which had preferred models that embody the early phase of clustering, with 
linearized matter fluctuations - since it is this final, nonlinear stage of clustering during 
which virialization occurs via the generation of vorticity and velocity dispersion, and hence 
represents the more astrophysically reasonable source of substantial causal backreaction. 

Noting that the only problem still remaining for this new concordance is the somewhat 
excessively large amount of clustering required to achieve it - that is, — 4, rather than 
what we consider to be more reasonable values like ~ 2 — 3 - we then determined that 
this problem could be successfully fixed (i.e., a good concordance generated with < 3) 
by introducing "early saturation", in which the clumping evolution function \l/(t) reaches its 
ultimate value of ^ somewhat in the past (z ~ 0.25), and then changes little thereafter. This 
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is a highly reasonable adjustment to the formalism, since in the real universe "gastrophysics" 
feedback exists which creates superheated baryons, sending large amounts of material back 
into the intergalactic medium, thereby slowing down the continued clustering of matter at 
late times; not to mention the likely slowdown of clustering due to the feedback effects of 
the backreaction, itself. 

The only major drawback of this new feature is the addition of an extra model input 
parameter - the epoch of saturation, £ Sat - which results in a degeneracy within (zsat> ^o) _ 
space, providing a range of models that all fit the Type la supernova data well, yet lead 
to significant differences for certain output cosmological parameters. The greatest variation 
in the output results due to this degeneracy occurs for the observable jerk parameter, j° bs , 
hence implying a loss of predictability for j° bs by our c ausal backreact ion formalism. This 



is a significant loss, given the previous findings from iBochnerl ( 1201 ll ) (without recursive 
nonlinearities) which had indicated that j® hs » 1 was the most distinctive signature of 
causal backreaction, thus serving as the clearest way for distinguishing it from Cosmological 
Constant ACDM (or from anything close to it), since flat ACDM always requires jo = 1. It 
thus becomes more difficult to find a falsifiable test of the causal backreaction paradigm, a 
test that is needed to definitively distinguish it from Dark Energy in order to eventually rule 
out one cosmological approach in favor of the other. 

Finally, concerning the 'ultimate' fate of the universe, we note that the incorporation 
of recursive nonlinearities tends to shut down any strong apparent acceleration effects fairly 
quickly once the ongoing clustering (i.e., the continued growth of finally stops. Even 

more dramatic is the way in which the metric perturbation function, J RNL (t), becomes essen- 
tially locked in place when approaching too close to unity, making it an even greater obstacle 
in terms of preventing the acceleration (apparent or otherwise) from completely taking over 
the cosmic evolution. This makes the scenario of a p erpetual, 'eternal' acceleration seem 



even less likely than it already did in IBochnerl ( 120111 ); though the now-unbounded nature 



of ^(t) could potentially provide some aid in producing a long-term acceleration, as long as 
virialized structure can continue to form on ever-larger cosmic scales, without any fundamen- 
tal upper limit to the sizes of coherent structures. Furthermore, the question of the ultimate 
cosmic fate is once again complicated by the possible backreaction contributions of gravita- 
tionally nonlinear terms, and the (unavoidable) eventual breakdown of the approximation 
of the universe as "smoothly-inhomogeneous" - both complications representing scenarios 
which our toy-model formalism is not presently designed to account for. 

In summary, we conclude that our causal backreaction formalism remains successful at 
generating an alternative cosmic concordance for a matter-only universe, without requiring 
any form of Dark Energy; though the necessary incorporation of recursive nonlinearities into 
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these models implies that a significantly stronger amount of such backreaction than before is 
now needed, acting throughout the crucial 'acceleration epoch' of z ~ 0.2 — 2 or so, in order 
to provide a degree of observed acceleration sufficient to match the cosmological standard 
candle observations. 



I am grateful to William Chan for computational support. 
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(r,0 = (()./„) 



Fig. 1. — Geometry for computing the inhomogeneity-perturbed metric at each point along 
the integrated path of a light ray from a supernova to our observation point at Earth. 




Fig. 2. — Residual Hubble diagrams for the causal backreaction model ^un{t) with 
(zinit, ^o) = (25,1.0), where the upper solid line represents the old version of the simu- 
lated cosmology without recursive nonlinearities (RNL), and the lower solid line represents 
the new version with RNL. Shown along with them for comparison (broken lines) are the 
flat SCDM and (Unionl-best-fit, = 0.713) Concordance ACDM cosmologies. 



-46- 




Fig. 3. — Residual Hubble diagrams for the two 'best' < 1 runs (as described in the 
text), selected from the new simulations with RNL but with the choice of model input 
parameters restricted to those used in BBI; specifically, these $md and ^ Sqr curves both have 
(zinit, ^o) = (5, 1.0). Also plotted here are the Unionl-best-fit flat SCDM and Concordance 
ACDM (fi A = 0.713) cosmologies. Shown along with these curves are the SCP Unionl SNIa 
data, here binned and averaged for visual clarity (bin size ALog 10 [l + z] — 0.01). Each 
theoretical model is individually optimized in H® hs to minimize its Xfh with respect to the 
Unionl SNIa data set; and for simplicity, instead of moving the SNIa data up or down 
for each different optimized if° bs value, the optimization is depicted here by plotting the 
residual Hubble diagram of the SNIa data versus a coasting universe of a single, fixed Hubble 
constant (i?° bs = 72 km s _1 Mpc _1 ), and then displacing each theoretical curve vertically, 
relative to the SNIa data, as appropriate for each fit. 
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Fig. 4. — Residual Hubble diagrams for three of the 'best' runs with > 1 (as described 
in the text), using model input parameters (solid lines increasing from lowest to highest): 
\l/s qr with (zinit, V&o) = (25,4.1) (curve almost indistinguishable from ACDM); \& L in with 
(zi n it, ^o) — (3,3.3); and $md with (z hlit , \& ) — (2,2.9). Shown along with them for com- 
parison (broken lines) are the flat SCDM and (Unionl-best-fit, Q\ — 0.713) Concordance 
ACDM cosmologies. 
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Fig. 5. — Residual Hubble diagrams for the early-saturation causal backreaction model 
^Sat(^) with (zinit, ^o) = (25,4.1), where the solid line peaking at z° hs ~ 1 — 2 (the curve 
almost indistinguishable from ACDM) represents a simulated cosmology with RNL but re- 
stricted to ^sat = (thus rendering it equivalent to ^sqr); and where the solid line peaking 
at z° hs ~ 10 — 20 represents the shift to z Sat = 3. Shown along with them for comparison 
(broken lines) are the flat SCDM and (Unionl-best-fit, fi A = 0.713) Concordance ACDM 
cosmologies. 
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Fig. 6. — Residual Hubble diag rams for early-saturation ^sat^) models with fixed zi n a — 25 
and varying zg at , with \l/o re-optimized for each chosen zgat value. From lowest to highest, 
the solid lines represent the runs with model input parameters: (zsat, ^o) — (0,4.1) (curve 
almost indistinguishable from ACDM); (z Sat , * ) = (0.25,2.6); (^sat,*o) = (0.5,2.3); and 
(^ S at,^o) = (1.0,2.2). Also plotted here (broken lines) are the Unionl-best-fit flat SCDM 
and Concordance ACDM (Qa = 0.713) cosmologies. Shown along with these curves are 
the binned and averaged SCP Unionl SNIa data (bin size ALog 10 [l + z\ — 0.01). Each 
theoretical curve is displaced vertically, relative to the SNIa data, to depict its individualized 
iJ° bs -opt imizat ion . 
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Fig. 7. — The clumping evolution functions are compared for several different early- 
saturation models, all with z Init = 25, plotted versus z FRW = [(t FRW /t FRW ) 2/3 - 1]. In 
order of the positions of the functions' 'corners' from left to right (increasing zs a t), the lines 
depict ^sat functions with the parameters: (zs a t, ^o) = (0,4.1); (^sat^o) = (0.25,2.6); 
(*Sat,*o) = (0.5,2.3); and (z Sa t,*o) = (1.0,2.2). 
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Table 1: Supernova Fit Quality and New- Concordance Success for Optimized RNL Runs 



z Init 


^o,o P t a 


2 b 
AFit 


Avg. % Dev. c 


$md Clumping Model Runs 


1 


2.8 


313.0 


7.5 


1.5 


2.8 


314.7 


4.0 


2 


2.9 


315.1 


5.8 


3 


3.1 


315.5 


14.1 


5 


3.6 


317.0 


30.4 


10 


4.3 


328.5 


52.2 


15 


4.1 


339.5 


51.5 


20 


4.1 


345.8 


52.6 


25 


4.1 


349.5 


53.1 


^Lin Clumping Model Runs 


1 


3.2 


311.0 


7.8 


1.5 


3.1 


312.8 


5.6 


2 


3.2 


313.6 


2.9 


3 


3.3 


314.0 


4.8 


5 


3.6 


314.3 


11.9 


10 


4.0 


314.7 


21.4 


15 


4.3 


314.9 


26.7 


20 


4.5 


315.0 


30.1 


25 


4.7 


315.1 


33.0 


^Sqr Clumping Model Runs 


1 


4.4 


311.5 


15.6 


1.5 


4.0 


310.4 


12.0 


2 


3.9 


310.5 


10.0 


3 


3.9 


311.0 


7.6 


5 


4.0 


311.4 


5.4 


10 


4.1 


311.7 


3.4 


15 


4.1 


311.8 


2.9 


20 


4.1 


311.8 


2.6 


25 


4.1 


311.8 


2.5 



a Each ^Opt chosen to minimize xlit f° r a gi ven ^init and \I>(i) function. 
b xlit computed versus the SCP Unionl SNIa data set. 

c Average Percent Deviation from the "Cosmic Concordance", defined (as described in the text) via cosmo- 
logical parameters obtained from the Unionl-best-fit ACDM model = 0.713 = 1 — Om, Xfh = 311.9). 
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Table 2: Output Cosmological Parameters from our 'Best' Runs with Recursive Nonlinearities 



2 Init 




Xlit a 


P b 

" Fit 


7 c 2° bs d H Obs e i/ FRW f £° bs 8 


QFRW h 


„,.Obs 
w 


„Obs 

Jo 


/Obs 
L A 


$md Clumping Model Runs 


1.5 


2.8 


314.7 


0.324 


0.53 1.16 69.54 41.96 13.95 


0.948 


-0.639 


0.59 


297.8 


2 


2.9 


315.1 


0.319 


0.62 1.17 69.58 38.62 14.44 


1.162 


-0.639 


0.51 


288.7 


xE'Lin Clumping Model Runs 


2 


3.2 


313.6 


0.340 


0.53 1.16 69.75 42.10 13.88 


0.946 


-0.675 


0.89 


297.3 


3 


3.3 


314.0 


0.334 


0.62 1.16 69.68 38.72 14.30 


1.159 


-0.667 


0.83 


287.2 


^Sqr Clumping Model Run 


25 


4.1 


311.8 


0.367 


0.53 1.14 70.07 42.32 13.64 


0.943 


-0.751 


1.73 


294.5 




Comparis 


on Values from Best-Fit 1 flat ACDM Model (O a 


= 0.713 


= i-n 


m) 








311.9 


0.380 


1.0 69.96 69.96 13.64 


0.287 


-0.713 


1.0 


285.4 


Comparison Values from Best-Fit •> _/?ai SCDM Model (fl\ = 0, 17 M = 1) 






608.2 


3.4E-22 


1.0 61.35 61.35 10.62 


1.0 


0.0 


1.0 


287.3 



a Xp it computed versus the SCP Unionl SNIa data set. 

b Each likelihood probability Ppit is derived from the corresponding Xfh using the Xn d f distribution with 
NuoF degrees of freedom, where AdoF = 304 for our $md, *Lin, and \&Sqr clumping models, Ad f = 305 
for the flat ACDM model, and AW = 306 for flat SCDM. 

c The integrated (Newtonian) gravitational perturbation potential at to, as modified by RNL, computed via 
Equations H51H31 
d Each z obs quoted here corresponds to z FRW = 1. 

e The HQ hs value (given here in km s _1 Mpc _1 ) for each run is found by minimizing its Xmt w ith respect to 
the SCP Unionl SNIa data set. 
•^Each Hq^ w is computed relative to the corresponding optimized HQ hs value for that run. 
9 All t$ hs values are listed here in GYr, and computed assuming no radiation (i.e., Qr = 0). 
h Al\ ri™ W va hies given here for the $md, ^Lin, and ^sqr models are normalized to ^S bs = 0.27. 
1 "Best-Fit" for the flat ACDM Model refers here to an optimization over and HQ hs . 
J "Best-Fit" for the flat SCDM Model refers here to an optimization over HQ hs . 
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Table 3: Output Cosmological Parameters from our RNL Runs with 'Early Saturation' 



^Sat 


*o,o P t 8 


2 

AFit 


P Fi t b 


T ^Obs frObs f/FRW 




qFRW 


^Obs 


•Ob s 

Jo 


;Obs 
( A 










^Sat Clumping Model Runs, 


Zlnit = 25 













4.1 


311.8 


0.351 


0.53 1.14 70.07 42.32 


13.64 


0.943 


-0.751 


1.73 


294.5 


0.25 


2.6 


313.5 


0.326 


0.58 1.15 69.60 40.24 


14.00 


1.054 


-0.620 


0.15 


289.7 


0.5 


2.3 


316.6 


0.284 


0.68 1.15 69.40 36.32 


14.65 


1.338 


-0.585 


-0.14 


279.8 


1 


2.2 


320.2 


0.238 


0.80 1.14 68.77 29.54 


15.75 


2.086 


-0.488 


-0.94 


259.9 



Comparison Values from Best-Fit flat ACDM Model (Qa = 0.713 = 1 - fl M ) 



311.9 0.380 ••• 1.0 69.96 69.96 13.64 0.287 -0.713 1.0 285.4 



a Each ^o,Opt is chosen to minimize xlit f° r a gi ven ^Sat, with zi^t — 25 and vE'(i) = ^Sat(i)- (Optimizations 
are performed with respect to the SCP Unionl SNIa data set.) 

b Each likelihood probability Ppit is derived from the corresponding xlit using the Xat DoF distribution with 
A^DoF degrees of freedom, where the inclusion of optimizable parameter zs a t now gives A^doF = 303 for our 
^Sat clumping models (even for zsat = 0), while A^oF = 305 is still used here for the flat ACDM model. 



